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1.  SUMMARY 


In  the  aviation  community,  there  is  a  high  priority  to  develop  all-electric  aircraft.  Electro¬ 
mechanical  actuation  (EMA)  systems  would  replace  traditional,  large,  heavy  and  difficult-to- 
maintain  hydraulic  actuation  systems.  This  movement  from  hydraulic  actuation  to  electrical 
actuation  enhances  the  flexibility  to  integrate  redundancy  and  emergency  system  in  future 
military  aircraft.  Elimination  of  the  hydraulic  fluid  removes  the  possibility  of  leakage  of 
corrosive  hydraulic  fluid  and  the  associated  fire  hazard,  as  well  as  environmental  concerns.  The 
switch  from  hydraulic  to  electrical  actuation  provides  additional  benefits  in  reduced  aircraft 
weight,  improved  survivability  and  improved  maintainability. 

Given  the  need  for  predicting  dynamically  generated  heat  and  power  demand,  an  accurate 
simulation  model  is  very  valuable  in  designing  an  EMA  system.  Advanced  analysis  techniques 
such  as  finite-element  method  (FEM)  and  intelligent  control  systems  such  as  field-oriented 
control  (FOC)  are  used  to  better  understand  the  sources  of  the  heat.  EMA  systems  are  complex 
and  testing  them  is  a  very  involved  procedure.  Such  systems  are  also  expensive  which 
discourages  testing  them  to  their  limits,  and  many  different  scenarios  are  needed  to  properly  test 
an  EMA.  Furthermore,  there  are  significant  details  which  cannot  be  practically  measured 
experimentally.  Therefore,  it  is  advantageous  to  simulate  the  EMA  with  only  its  key  parameters. 
The  EMA  simulation  model  presented  in  this  report  meets  this  goal. 

We  address  the  two  most  important  components  of  the  EMA  -  electric  motor  and  power  drive. 
The  model  takes  a  flight  mission  (a  desired  position  profile  and  a  load  force  profile)  and  an 
environmental  temperature  profile  as  input  and  predicts  the  thermal  behavior  of  the  electric 
motor  and  the  transient  power  demands  of  the  EMA.  It  includes  a  position  field-oriented  control 
(FOC),  a  pulse-width  modulation  (PWM)  component,  an  electro-mechanical  dynamics 
component,  and  a  thermal  component.  Model  parameters  are  obtained  by  experimentation  or 
FEM  models.  The  key  parameters,  including  nonlinear  inductance,  rotor  flux  linkage,  and 
thermal  resistances  and  capacitances,  are  tuned  using  FEM  models  of  a  commercial  actuator. 

The  full  simulation  model  is  a  nonlinear,  lumped-element  model  (NL-LEM).  The  model 
corresponds  to  a  linear  actuator  with  a  permanent  magnet  (PM)  rotational  motor  connected  to  a 
rotational-to-linear  drive  train.  It  was  shown  that  the  heat  generation  profiles  are  indeed  highly 
dynamic,  and  control  plays  a  powerful  role  in  how  much  excess  heat  is  generated  and  how  much 
the  power  demand  spikes. 

For  the  thermal  component,  a  comprehensive  FEM  model  has  been  developed  to  predict  the 
temperature  distribution  within  the  electric  motor  using  the  heat  load  generated  in  a  given  flight 
profile.  However,  the  FEM  model  needs  tremendous  computational  resources.  A  reduced-order 
thermal  model  (lumped-node  thermal  network)  is  validated  and  used  instead,  which  can  reduce 
computational  resources  considerably.  The  network  can  simulate  the  temperature  time  histories 
of  the  various  components  in  the  EMA  electric  motor  during  a  flight  mission. 

The  thermal  lumped-node  model  is  combined  with  the  electromagnetics  portion  mentioned 
above  to  dynamically  calculate  the  temperature  within  the  electric  motor.  FEM  results  obtained 
with  or  without  radiation  are  compared  with  the  results  from  the  lumped-node  model.  Radiation 
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is  important  when  convection  cooling  is  insufficient  and  the  motor  temperature  is  high.  In 
situations  when  the  heat  is  excessive  and  cannot  be  removed  from  the  motor  real-time,  we  also 
demonstrate  that  phase  change  materials  (PCM)  is  very  effective  in  keeping  the  motor 
temperature  near  the  melting  point  of  the  PCM  until  all  PCM  has  melted. 

We  also  study  the  thermal  problem  of  the  control  and  drive  units  of  EMAs,  and  build  a  model  to 
calculate  and  simulate  the  power  loss  and  heat  generation  in  the  driver  board.  The  driver  unit 
consists  of  a  power  inverter,  a  power  dissipating  resistor  and  a  control  circuit.  The  power  loss  of 
each  part  is  studied.  The  heat  loss  in  the  power  inverter  comes  mainly  from  the  power  switches  - 
IGBTs.  The  on-state  loss  is  proportional  to  the  current  of  the  motor,  and  the  switching  loss  is 
determined  by  the  switching  frequency  as  well  as  currents  and  voltages.  The  power  loss  in  the 
power  dissipating  resistor  is  determined  by  the  regenerative  power,  the  capacitor  and  the  control 
algorithm  to  stabilize  the  bus  voltage,  which  varies  from  different  mission  profiles  and  different 
applications.  The  power  loss  on  the  control  circuit  is  negligible  compared  with  the  power  loss  on 
IGBTs  and  the  power  dissipating  resistor,  and  generates  very  little  heat. 

As  a  brief  summary,  this  report  presents  a  dynamic  heat  generation  model  and  a  thermal  model 
for  a  PM  motor  and  power  drive.  The  input  to  the  model  includes  a  flight  mission  profile  and  an 
environmental  temperature  profile,  along  with  the  key  parameters  unique  to  the  EMA.  The 
model  compares  well  with  the  detailed  finite  element  analysis.  The  model  was  also  used  to 
compare  its  temperature  predictions  with  temperature  measurements  of  the  motor  and  power 
drive  provided  by  the  Air  Force  Research  Laboratory.  Though  the  agreement  was  good,  due  to 
the  low  heat  load  associated  with  the  test  cases,  the  model  could  not  be  validated  conclusively. 
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2.  INTRODUCTION 


The  development  of  all-electric  aircraft  is  a  high  priority  in  the  avionics  community  [1].  Current 
aircraft  use  a  combination  of  hydraulic,  pneumatic,  and  electric  systems.  However,  future 
airplanes  are  expected  to  use  a  single,  electric  system,  with  electromechanical  actuators  (EMAs) 
replacing  hydraulic  pistons.  Such  a  system  would  reduce  the  cost  to  build,  operate,  and  maintain 
aircraft  [2],  It  would  also  make  aircraft  lighter,  more  reliable,  safer,  and  more  easily 
reconfigurable,  improving  the  turnaround  for  new  technology  [3], 

There  are  two  major  obstacles  in  replacing  hydraulic  actuators  with  EMAs:  heat  generation 
management  and  power  management.  Firstly,  unlike  hydraulic  actuators,  EMAs  do  not  have  the 
inherent  advantage  of  recirculating  fluid  to  cool  their  components  [4],  Rather,  windings  in  EMAs 
can  overheat  rapidly  depending  on  the  demands,  at  which  point  materials  can  degrade. 
Neodymium-iron-boron  is  one  of  the  most  powerful  permanent  magnets  for  electric  motors,  but 
it  begins  to  demagnetize  at  relatively  low  temperatures  [5],  with  an  operating  range  of  120°C  to 
180°C.  Although  effective,  liquid  cooling  reinstitutes  one  of  the  systems  that  are  being 
eliminated.  Primary  flight  control  actuators  are  of  particular  relevance  because  they  are 
continually  engaged  during  flight,  are  required  to  accelerate  rapidly,  and  are  often  faced  with 
high  wind  loads.  As  the  complexities  of  aircraft  systems  exponentially  grow  and  the  demand  for 
lightweight  composites  increases,  thermal  concerns  are  increasingly  significant.  Heat  generation 
in  motor  windings  can  be  highly  dynamic  and  localized.  Predicting  how  much  heat  would  be 
generated  by  an  EMA  is  critical. 

Secondly,  EMAs  have  the  advantage  of  only  using  power  when  active,  whereas  hydraulic 
systems  require  that  a  pump  continuously  maintain  line  pressure.  This  should  help  to  make  EMA 
systems  more  efficient.  However,  this  advantage  brings  with  it  a  major  complication:  the  power 
draw  from  the  electrical  supply  is  also  highly  dynamic.  Compounding  this  problem  are  the 
highly  variable  demands  of  all  the  electrical  systems  of  an  entire  aircraft  drawing  power  from 
one  source.  Predicting  how  much  power  will  be  drawn  during  operation  is  critical  for  sizing 
system  components  and  choosing  the  power  supply  type  (e.g.  battery-super-capacitor  hybrid). 

In  the  next  two  chapters,  we  will  address  the  heat  generation  and  power  management  of  the  two 
most  important  components  of  the  EMA  -  electric  motor  and  power  drive.  The  model  takes  a 
flight  mission  (a  desired  position  profile  and  a  load  force  profile)  and  an  environmental 
temperature  profile  as  input  and  predicts  the  thermal  behavior  of  the  electric  motor  and  the 
transient  power  demands  of  the  EMA. 
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3.  DYNAMIC  HEAT  GENERATION  MODELING  OF  HIGH 
PERFORMANCE  ELECTROMECHANICAL  ACTUATOR 

3.1  Introduction 

Given  the  need  for  predicting  dynamically  generated  heat  and  power  demand,  an  accurate 
simulation  model  is  tremendously  valuable  in  designing  an  EMA  system.  EMA  systems  are 
complex  and  testing  them  is  an  involved  procedure.  Such  systems  are  also  expensive  which 
discourages  testing  them  to  their  limits,  and  many  different  scenarios  are  needed  to  properly  test 
an  EMA.  Furthermore,  there  are  significant  details  which  cannot  be  practically  measured 
experimentally.  Therefore,  it  is  advantageous  to  simulate  the  EMA  with  only  its  key  parameters. 

The  EMA  simulation  model  presented  in  this  chapter  meets  this  goal.  It  includes  a  position  field- 
oriented  control  (FOC),  a  pulse-width  modulation  (PWM)  component,  an  electro-mechanical 
dynamics  component,  and  a  thermal  component.  Model  parameters  are  obtained  by 
experimentation  or  finite  element  method  (FEM).  The  full  simulation  model  is  a  nonlinear, 
lumped-element  model  (NL-LEM).  The  model  corresponds  to  a  linear  actuator  with  a  permanent 
magnet  (PM)  rotational  motor  connected  to  a  rotational-to-linear  drive  train  (e.g.  a  ball  screw  or 
a  roller  screw). 

The  methodology  presented  here  provides  a  means  of  predicting  heat  generation  and  the  overall 
power  demands  of  an  EMA.  Such  tools  will  help  in  reaching  the  goal  of  an  all-electric  aircraft. 

Simulation  results,  including  heat  generation  and  power  demand,  for  a  commercially  available, 
test  EMA  are  shown  in  the  latter  part  of  this  chapter. 

3.2  Model  Layout 

The  four  main  components  of  the  NL-LEM  are  shown  in  Figure  1. 
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Figure  1  -  Model  Layout  Showing  the  Connection  of  the  four  Major  Components:  Field- 
Oriented  Control,  Pulse-Width  Modulation,  Electromechanical  Dynamics,  and  the 

Thermal  Component 

The  model  takes  the  EMA  rod’s  desired  motion  profile,  x* ,  as  an  input  to  the  field-oriented 
control  (FOC).  Based  on  the  present  state  of  the  direct  and  quadrature  (DQO)  currents,  id,iq ,  and 

the  actual  position  of  the  EMA  rod,  x ,  the  FOC  algorithm  calculates  the  reference  voltages, 
u/,u * .  The  PWM  component  takes  these  reference  voltages  as  well  as  the  present  state  of  the 
DC  supply  voltage,  uDC ,  and  calculates  the  PWM  voltages,  ud,uq .  These  are  the  actual  input 
voltages  to  the  electromechanical  dynamics.  With  these  and  the  present  load  force,  FL ,  the 
electromechanical  dynamical  equations  lead  to  the  calculation  of  the  currents,  id,iq ,  EMA 
position,  x,  and  dynamically  generated  heat,  Q.  An  output  of  the  electromechanical  component  is 
the  DC  current  to  the  EMA  system,  iDC .  Together  with  the  DC  bus  voltage,  it  represents  the 
power  demands  of  the  EMA  system  at  its  terminals.  The  heat  generated  in  the  motor  is  the  input 
to  the  thermal  component.  With  the  present  value  of  the  environmental  temperature,  Tenv ,  the 
thermal  equations  lead  to  the  nodal  temperatures,  7] ,  and  the  updated  value  of  resistance,  Rs , 
which  is  fed  back  into  the  electromechanical  dynamics.  The  nodal  temperatures  represent  the 
thermal  behavior  of  the  EMA  system.  Table  1  summarizes  many  of  these  terms. 
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Table  1  -  Simulation  States 


Symbol 

Name 

Units 

X 

Actual  position  of  EMA  rod 

mm 

X 

Desired  position  of  EMA  rod 

mm 

dme 

Actual  rotor  angle  scaled  by  P/2 

rad 

oj 

Desired  rotor  angle  scaled  by  PI 2 

rad 

V 

Actual  velocity  of  EMA  rod 

mm/s 

V 

Desired  velocity  of  EMA  rod 

mm/s 

fl 

Load  force 

N 

Fm 

Motor-generated  force 

N 

Fg 

Gravitational  force 

N 

Ff 

Frictional  force 

N 

ud 

Actual  direct  voltage 

V 

u 

q 

Actual  quadrature  voltage 

V 

«/ 

Desired  direct  voltage 

V 

u* 

Desired  quadrature  voltage 

V 

Actual  direct  current 

A 

Actual  quadrature  current 

A 

id 

Desired  direct  current 

A 

c 

Desired  quadrature  current 

A 

UDC 

DC  bus  voltage 

V 

lDC 

DC  bus  current 

A 

Q 

Generated  heat 

W 

Rs 

Phase  resistance 

Q 

Tenv 

Environmental  temperature 

°C 

Ti 

Nodal  temperature 

°C 

Notice  that  the  desired  position,  x* ,  the  load  force,  FL ,  the  environmental  temperature,  Tenv ,  and 
the  DC  voltage,  uDC ,  all  can  be  variable  inputs  to  the  model.  An  EM  A  can  be  run  through  a 
motion  profile  in  the  presence  of  a  variable  load  profile  while  the  electrical  power  demand  and 
the  temperatures  are  calculated. 

3.3  Parameter  Measurement  and  FEM  Calculation 

Table  2  shows  a  summary  of  the  key  EMA  parameters.  Some  of  the  parameters,  such  as  the 
number  of  slots  and  number  of  poles,  are  easily  obtained  from  the  actuator’s  documentation,  but 
others  must  be  measured  or  carefully  calculated.  Some  of  the  measurements  are  performed  while 
the  actuator  is  coupled  with  a  hydraulic  press  or  some  other  active  load. 
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Table  2  -  EMA  Parameters 


Symbol 

Name 

Units 

K 

Gearing  ratio 

rad/m 

s 

Number  of  slots 

p 

Number  of  poles 

Ff 

f,  max 

Maximum  power  train  friction 

N 

m 

Rod  mass 

kg 

I 

Rotor  moment  of  inertia 

kg-mm2/rad2 

R s,ref 

Phase  resistance  at  t  , 

ref 

Q 

Tref 

Reference  temperature 

°C 

aR 

Temperature  coefficient  of  resistance 

1/°C 

^ PM 

PM  flux  linkage 

Wb 

aMJ 

Direct  inductance 

mH 

A  (A  a) 

Quadrature  inductance 

mH 

3.3.1  Measurement 

a.  Gearing  Ratio,  Ncr 

The  gearing  ratio  is  carefully  calculated  by  counting  the  teeth  on  each  gear  and  measuring  the 
linear  displacement  of  the  actuator  rod  for  one  full  revolution  of  the  rotor. 

b.  Power  Train  Friction,  Ff 

The  hydraulic  actuator  is  coupled  to  the  unpowered  EMA  and  drives  the  EMA’s  rod  up  and 
down.  The  static  friction  of  the  EMA’s  drive  train  is  obtained  by  using  the  force  and  stroke  data 
collected  during  this  test,  and  the  mean  forces  required  to  push  or  pull  the  EMA  rod  are 
calculated.  The  magnitudes  of  these  forces  are  unequal  due  to  a  bias.  This  bias  exists  because  of 
the  weight  of  the  rod  itself,  and  the  difference  between  this  bias  and  the  push  and  pull  forces  is 
the  friction  force. 

c.  Rod  Mass,  m 

The  mass  of  the  EMA’s  entire  rod  is  derived  from  the  friction  test.  The  bias  in  the  force  required 
to  move  the  EMA  rod  is  due  to  the  weight  of  the  rod.  This  translates  to  a  certain  mass. 

d.  Rotor  Moment  of  Inertia,  I 

The  rotor  moment  of  inertia  is  obtained  by  measuring  the  dimensions  of  the  rotor  and  the 
material  densities  and  calculating  the  moment  of  inertia. 

e.  Phase  Resistance,  Rs 

The  phase  resistance  is  measured  by  connecting  leads  to  two  of  the  three  phase  lines  into  the 
motor  and  measuring  the  open-circuit  line-to-line  resistance.  This  value  is  then  divided  by  two  if 
the  motor  is  a  wye-connected  motor. 
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3.3.2  FEM  Calculation 


The  material  composition  and  geometry  of  the  motor  are  used  in  designing  an  FEM  model  in 
Ansoft’s  Maxwell  2D.  Values  of  rotor  angle,  torque  angle,  and  current  amplitude  are  varied  as 
multiple  simulation  runs  are  performed. 

Many  of  the  motor  modeling  researchers  in  the  past  [6],  [7]  and  even  present  [8],  [9]  use  constant 
parameter  values  in  simulations,  although  some  use  nonlinear  inductances  [10],  [11],  [12],  For 
this  chapter,  FEM  is  used  to  obtain  Ld  and  Lq  as  nonlinear  functions  of  id  and  iq .  The 

derivatives  of  inductance  with  respect  to  current  are  numerically  calculated  from  these  tables. 

The  inductances  and  their  derivatives  with  respect  to  DQO  currents  are  stored  in  tables  which  are 
then  indexed  during  simulation. 

The  FEM  model  is  simulated  with  currents  well  above  the  motor’s  normal  operating  range  in 
order  to  capture  the  behavior  during  highly  transient  moments.  The  torque  angle  should  be  varied 
from  0  to  7i  radians  to  capture  the  full  character  of  the  inductances.  Even  if  the  rotor  of  an  EMA 
is  round,  the  rotor  angle  has  an  effect  on  the  DQO  reference  frame  inductances.  Therefore,  the 
rotor  angle  should  also  be  varied  over  one  period  of  variation  in  the  DQO  inductances.  Such 
variation  occurs  mostly  at  relatively  high  current  values  (i.e.  at  saturation),  so  this  variation  is 
averaged  out,  but  the  calculation  of  inductance  at  a  single  rotor  angle  cannot  be  assumed  to  be 
representative. 

Note  that  the  peak  inductance  values  do  not  necessarily  occur  at  zero  current,  nor  are  the 
inductances  mirrored  about  the  quadrature  axis.  An  example  of  this  is  in  Figure  2.  This  should 
clarify  the  highly  variable  nature  of  inductances  and  the  need  to  use  nonlinear  modeling  instead 
of  assuming  constant  values. 


Figure  2  -  Direct- Axis  Inductance  Along  the  id  Axis  (blue)  and  the  Inductance  Value  at 
Zero  Current  (Green  Circle).  The  Inductance  is  not  mirrored  about  the  Quadrature  Axis 
(going  into  the  page),  so  its  Values  for  Positive  id  do  not  match  those  for  Negative  id 

3.4  Electromechanical  Dynamics 

The  motor  dynamics  are  modeled  by  four  primary  dynamical  equations  in  the  DQO  reference 
frame  [13].  The  first  two  of  these  equations  are 
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did  _Lqfid~Ldfi± 


(1) 


dt  d,ddLqq  LdqLqd 


dig  _  LdJdq  LqdG^ 
dt  d,ddLqq  —  LdqLqd 


These  come  from  the  voltage  Eq.  3  and  Eq.4: 


di ,  di  dO 

u,  =  R  i,  +  L ..  -d-  +  L .  — — -Z,  z  — ^ 

'  "A  *  *  ’  5  dt 


u  =  R  i  +  L  .  +  L  — — 

"  s  ?  ^  *  w  dt 


+  (A/A  +  ) 


ddme 

dt 


(2) 

(3) 

(4) 


where 


^  -L  +i 

^ dd  ~  +  ld 

dij 


(5) 


_  dLj_ 

dq  d  di 


Ld=i^ 
““  “  dL 


dL 

L  =L  +i  — q- 

qq  q  q  di 


Gd=ud-Rsid+{Lqiq)^ 


G,  ~  11  u  RJg  (AT/  +  ^PM  ) 


Ome={^cr)- 


ddmt 

dt 


(6) 

(7) 

(8) 

(9) 

(10) 

(11) 


These  equations  carry  the  assumption  that  the  variations  in  inductances  over  rotor  angle  can  be 
reasonably  averaged  out,  that  the  flux  linkage  from  the  permanent  magnets  appears  only  in  the 
direct  axis,  and  that  the  motor  is  balanced. 

The  voltages  are  provided  by  the  PWM  component;  the  current  rates  are  calculated  using  Eq.  1 
and  Eq.  2;  and  the  currents  are  found  by  integration. 

In  [14],  there  is  a  detailed  analysis  of  the  mo  tor- generated  force  equation.  Given  the  assumptions 
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made  for  Eq.  1  -  Eq.  1 1 ,  the  motor-generated  force  equation  is 


Fm  ~NJq  Wm  +  (Ld  -  Lq  )d  ]  (12) 

The  third  dynamical  equation  is  the  force-speed  equation  (Eq.  13)  that  links  the  motor-generated 
force  to  the  acceleration  it  produces: 


dv 

dt 


+  FL+Fg+Ff) 


(13) 


where  J  is  the  combined  linear  inertia  from  all  moving  parts  in  the  EMA.  The  value  J  is 
calculated  from  the  mass  of  the  rod  and  the  angular  moment  of  inertia  by 


J  = 


I—N  2  +m 
2 


(14) 


As  long  as  velocity  is  zero,  the  frictional  force  is  defined  to  be  equal  and  opposite  to  the  sum  of 
the  other  component  forces  up  to  a  maximum  friction  value.  This  keeps  the  net  force  and, 
therefore,  velocity  at  zero.  However,  when  the  sum  of  the  other  component  forces  exceeds  this 
maximum  friction,  the  magnitude  of  the  friction  holds  at  this  maximum  value  and  its  sign  always 
opposes  the  direction  of  motion  Eq.  15. 


Ufm+fl+fs) 

■\fm+fl+fs\ 

|-sgn(v)-F/>max 

\fm+fl+f\ 

<  F, 

J  ,max 

>  F 

/,  max 


(15) 


Note  that  air  friction  is  not  included  because  EMAs  use  position  controlled  motors  which  do  not 
reach  excessive  speeds  where  air  friction  is  a  practical  consideration. 

Finally,  the  forth  dynamical  equation  is 


dx 

dt 


v 


This  is  simply  the  definition  of  velocity  and  is  integrated  to  get  position. 


(16) 


3.5  Field-oriented  Control 

A  thorough  analysis  of  dynamic  motor  heat  generation  should  include  motor  drive  dynamics. 
This  model’s  field-oriented  control  (FOC)  translates  the  desired  motion  profile  into  simulation 
drive  voltages.  The  control  is  an  adaptive  closed-loop  control,  which  is  a  type  of  feedback 
linearization  of  the  speed  and  position  errors. 

In  the  following  derivation,  the  star  superscript  is  used  to  indicate  desired  values.  Starting  with 
velocity  and  position  error,  a  desired  linear  acceleration,  a* ,  is  developed: 
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a  -kt 


(17) 


a 


where  A tc  is  the  time  step  of  the  control  loop  and  ka  and  kv  are  control  coefficients.  This 
equation  tends  toward  changing  the  present  velocity  and  position  to  match  the  desired  velocity 
and  position.  However,  this  should  happen  only  at  a  certain  rate  and  not  necessarily  all  in  one 
time  step. 

The  desired  velocity  must  be  redefined  based  on  this  newly  calculated  desired  acceleration: 


v  =v  +  a  A tc 


(18) 


If  the  load  force  is  predicted  by  estimating  the  load  force  from  the  previous  step  and  if  the 
frictional  and  gravitational  forces  are  neglected,  then  Eq.  13  becomes  Eq.  19. 


p  = 

-L  A/T 


aJ~Fg-Ff 


3  P 

NCr  ~Z  A PM 

Z  Z  y 


~Fg~Ff 


FL 


Solving  Eq.  12  for  i  and  setting  ij  to  zero  yields 


(19) 


3  P 

-Ncr-A,PM 

2  2 


Substituting  Eq.19  into  Eq.20  results  in  Eq.21. 


-  (a  -tty  u 


3  P 

-n„—; l 


The  currents  and  the  current  rates  are  then  limited. 
Lastly,  the  desired  voltages  are  calculated: 


*  dir  *  dO 

-R^  -v.  ~t 


(20) 


(21) 


(22) 


=  Rsiq 


+  L 


q  0 


di 

q 

dt 


+  +  ^PM 


yld,„t 

dt 


(23) 


where  Ld0  is  the  direct  inductance  at  zero  current,  Lq0  is  the  quadrature  inductance  at  zero 
current,  and  9m*  is  the  desired  mechanical  angle  multiplied  by  the  number  of  pole  pairs  PI  2 .  Its 
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relation  to  desired  position  is 


(24) 


3.6  PWM  Component 

The  PWM  component  of  the  simulation  model  uses  a  trailing  edge  modulation  carrier,  as 
demonstrated  in  Figure  3.  The  voltages  are  translated  into  the  ABC  reference  frame.  The  PWM 
algorithm  is  performed  on  those  voltages.  Then,  these  pulsed  voltages  are  translated  back  into  the 
DQO  reference  frame. 


Comparison  of  input  and  carrier 


Figure  3  -  The  General  Concept  of  the  PWM.  A  Trailing  Edge  Carrier  Wave  (top)  is 
compared  to  a  Normalized  Sinusoidal  Phase  Voltage.  When  the  Normalized  Wave  is 
Greater  than  the  Carrier  Wave,  the  Output  PWM  Voltage  is  turned  on  (Bottom). 

Otherwise,  it  is  off 


3.7  Inputs 
a.  Scaling 

The  stroke  and  load  force  data  are  recorded  from  actual  flight  profiles.  However,  these  profiles 
correspond  to  an  actuator  with  performance  characteristics  that  might  not  match  those  of  a  test 
actuator.  The  stroke  and  force  profiles  are  scaled  for  a  range  which  is  compatible  with  the  test 
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actuator.  The  stroke  profile  is  also  adjusted  to  match  the  acceleration  limits  of  the  actuator  by 
stretching  the  profile  over  time. 

b.  Force-Stroke  Equation 

In  the  absence  of  actual  flight  profiles,  a  force-stroke  relationship  can  be  useful  in  simulating  the 
force,  given  a  theoretical  stroke  profile.  This  force-stroke  equation  comes  from  analyzing  actual 
flight  data.  The  equation  is 


E  —  ip '  Cd  +Cd0)- 

[(.t  -  x0 )  +  C2  (x  -  x0  f  +  C,  (x  -  x0  )3  Js 2  + 

F0+cv  +  Ma  (25) 

q  _  j^dl  '  X  <  X0 
\Cd2  :x>x0 

where  p  is  the  air  density  which  can  be  derived  from  the  altitude  of  the  aircraft,  Cdl ,  Cd2 ,  Cd0 , 

C  2 ,  and  C  3  are  coefficients,  5  is  the  aircraft  speed,  F0  is  the  force  bias  due  to  the  weight  of  the 

aileron,  c  is  the  net  friction  coefficient,  Mis  the  effective  mass  of  the  aileron  and  EMA,  x  is  the 
stroke  position,  xo  is  a  stroke  bias,  v  is  the  first  derivative  of  x,  and  a  is  the  second  derivative  of 
x.  The  effect  of  air  density  depends  on  whether  the  flap  is  down  or  up  ( x  <  x0  or  x  >  x0 ).  The 
effect  of  aircraft  speed  is  dependent  on  a  Taylor  series  function  of  the  stroke. 

The  result  of  such  fitting  of  the  data  on  hand  is  shown  in  Figure  4  where  the  force  is  a  function  of 
stroke. 


Figure  4  -  Actual  Force  Profile  (Blue)  and  Calculated  Force  Profile  (Green) 

3.8  Thermal  Modeling  of  the  Motor 

High  performance  electric  actuation  systems  (HPEAS)  are  being  developed  for  new  generation 
flight  control  surfaces  under  the  Air  Force  Research  Laboratory’s  (AFRL)  Integrated  Vehicle 
and  Energy  Technology  (INVENT)  program.  The  new  distributed  electro-mechanical  actuation 
(EMA)  systems  would  replace  the  traditional,  large,  heavy  and  difficult-to-maintain  hydraulic 
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actuation  system.  The  electric-based  actuation  is  an  enabling  technology  to  simplify  the  whole 
power  generation  and  distribution  system. 

The  benefits  brought  from  this  transition  are  noticeable.  With  the  traditional  hydraulic  power 
system,  the  pumps  always  run  at  3000  psi,  whether  the  actuators  are  used  or  not.  All  the  energy 
into  the  hydraulic  system  becomes  heat  and  the  heat  dissipated  must  be  dumped  into  the  fuel 
system  and/or  discarded  to  the  ambient  air.  The  EMA  system  provides  power  on  demand.  When 
flying  level  at  near  constant  speed,  the  electrical  power  requirement  is  small.  Only  when  the 
aircraft  performs  turn,  pitch,  roll  and  yaw  in  high  acceleration  would  the  actuators  need  to 
provide  high  power  to  move  the  flight  control  surfaces.  The  new  electrical-based  actuation 
systems  generate  less  waste  heat  than  the  hydraulic  power  systems  [15]. 

3.8.1  Finite  Element  Analysis 

Figure  5  is  the  3-D  model  of  a  typical  Permanent  Magnet  Synchronous  Machine  (PMSM)  servo 
motor.  The  motor  design  features  are  implemented  in  a  3-D  frame  based  on  the  motor  structure 
with  a  potted  epoxy,  square  case  and  closed  no-venting  cover  plate.  The  motor  has  a  12-slot 
stator  and  a  10-pole  rotor  as  shown  in  Figure  6.  This  motor  is  designed  to  output  10  hp  electrical 
power,  with  an  efficiency  of  92%.  The  total  power  loss  in  the  motor  is  600  W. 


Figure  5  -  3-D  Model  of  the  PMSM  Motor  Design 
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Figure  6  -  Cross-Section  of  Motor  and  Rotor 

Due  to  the  symmetry,  the  geometry  of  the  model  is  divided  into  a  quarter  section  and  imported 
into  ANSYS  [16]  for  the  FEA  simulation.  Figure  7  shows  the  quarter  section  of  the  motor  where 
the  various  components  and  the  air  gap  are  shown.  The  total  thermal  load  is  given  in  Figure  8. 
The  estimated  distribution  of  the  power  loss  in  the  motor  components  is  listed  in  Table  3.  We 
assume  all  the  heat  loss  occurs  at  the  copper  winding  and  end  turns,  and  neglect  the  other  losses 
because  they  are  all  relatively  small.  In  the  computation  domain  the  power  loss  is  25%  of  total 
power  loss  600  W  or  150  W. 


1.  copper  winding 

3.  Iron  stator 

4.  case 

5.  case  surface 

6.  shaft 

7.  shaft  front  end 

8.  shaft  back  end 

9.  magnet 

10.  air  gap 

12.  front  bearing 

13.  rear  bearing 

14.  front  cover 


Figure  7  -  Computation  Domain 
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Figure  8  -  Step  Heat  Load  of  the  Motor 


Table  3-  Power  Losses  in  Motor  Components  in  a  Quarter  Section 


Power  loss  (W) 

Percentage  of  total  power  loss 

IC  (copper  loss) 

125.2 

83.5% 

ICF  (front  copper  end  turn  loss) 

8.1 

5.4% 

ICB  (rear  copper  end  turn  loss) 

8.1 

5.4% 

IS  (stator  loss) 

6.0 

4.0% 

IM  (magnet  loss) 

0.8 

0.6% 

IW  (windage  loss) 

0.4 

0.2% 

IBF  (front  bearing  loss) 

0.8 

0.6% 

IBB  (rear  bearing  loss) 

0.5 

0.6% 

ANSYS  generally  gives  temperature  information  in  the  whole  calculation  domain.  However,  in 
practice,  only  several  locations  are  needed  to  assess  the  thermal  behavior  of  the  motor.  We  focus 
on  the  temperatures  of  the  winding,  stator,  magnet  and  bearing.  It  is  noted  that  for  every 
component,  the  temperature  is  not  uniform.  The  temperature  range  can  be  characterized  by  the 
maximum  and  minimum  temperatures  within  the  particular  component.  For  the  copper  winding, 
these  two  temperatures  are  quite  close  because  of  the  high  thermal  conductivity  of  copper. 

While  in  the  stator  these  two  temperatures  differ  by  a  few  degree  Celsius.  An  external 
convective  heat  transfer  coefficient  of  h=100  W/m  K  (forced  convection),  and  heat  transfer 
coefficient  of  h=10  W/m  K  (natural  convection)  are  applied  to  the  case  of  the  motor.  Table  4 
shows  the  results  from  the  simulation  with  an  ambient  temperature  of  22°C.  The  node  number  in 
Table  4  corresponds  to  the  part  number  shown  in  Figure  7.  When  radiation  (with  a  surface 
emissivity  <r  of  0.85  assumed  for  the  motor  case)  is  applied  to  the  model  along  with  forced 
convection  (h=100  W/m2-K),  the  results  show  a  lower  temperature  (~20°C).  This  means  though 
forced  convection  is  the  dominant  mode  of  heat  transfer,  radiation  losses  should  be  included. 
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With  natural  convection  (h=10  W/m  K),  radiation  is  the  dominant  mode  of  heat  transfer,  but  the 
motor  temperature  is  much  too  high.  Thus  it  can  be  concluded  at  a  heat  load  of  600  W,  natural 
convection  is  not  an  acceptable  cooling  method.  Additional  discussions  on  radiation  modeling 
are  given  later  in  this  chapter. 


Table  4  -  Steady  State  Simulation  Results  with  AN SYS 


Node  name 

Material 

Node 

number 

Temp. 

(h=100)  in  °C 

Temp. 

(h=100)  in  °C 
w/  radiation, 
£=0.85 

Temp. 
(h=10)  in  °C 
w/  radiation, 
£=0.85 

Copper  winding 

Copper 

1 

162 

144 

396 

Iron  stator 

Iron 

3 

161 

143 

395 

Case 

Aluminum  alloy 

4 

152 

133 

385 

Case  surface 

5 

150 

130 

382 

Front  cover 

Aluminum  alloy 

14 

141 

124 

377 

Shaft 

Stainless  steel 

6 

149 

125 

384 

Shaft  front  end 

Stainless  steel 

7 

145 

124 

383 

Shaft  rear  end 

Stainless  steel 

8 

148 

124 

383 

Magnet 

NdFeB 

9 

149 

125 

385 

Air  gap 

Air 

10 

160 

135 

395 

Front  bearing 

Steel 

12 

143 

124 

381 

Rear  bearing 

Steel 

13 

147 

125 

381 

For  the  step  heat  load  shown  in  Figure  8,  the  steady-state  temperature  field  and  the  transient 
temperature  response  of  the  different  parts  inside  the  motor  are  given  in  Figures  9  and  10.  The 
initial  temperature  of  the  motor  is  uniform  at  22°C  and  ambient  temperature  is  also  at  22°C. 

The  heat  transfer  coefficient  is  100  W/m2K  and  radiation  is  included.  The  time  constants  for  the 
temperature  rise  of  the  motor  components  are  between  250-500  seconds. 
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G:  Steady-State  w/  Radiation,  h=100 

Temperature 
Type:  Temperature 
Unit:  -C 
Time:  1 


0.060  (m) 


Figure  9  -  Steady  State  Thermal  Simulation  Result  (Forced  Convection) 


Time  (s) 


Figure  10  -  Transient  Temperature  Response  of  the  Step  Heat  Load  Shown  in  Figure  8 

3.8.2  Lumped  Node  Model 

A  lumped  node  thermal  network  to  represent  the  temperature  of  every  solid  part  of  the  EMA  has 
been  developed.  The  thermal  resistances  and  capacitances  between  the  nodes  can  be  treated  as 
electrical  resistances  and  capacitances  (Table  5).  Hence  the  temperature  of  every  node  can  be 
solved  as  the  voltage  in  this  equivalent  network.  This  approach  is  well  developed  in  the  motor 
design  industry  [17,  18].  Some  commercial  software  packages  have  already  included  the  thermal 
network  simulation,  i.e.  Motor-CAD  [19].  Some  studies  have  been  made  with  FEA  analysis  and 
experimental  testing  has  shown  that  such  an  approach  is  valid  [20]. 
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In  the  lumped  node  network  we  use  a  single  temperature  to  represent  each  component,  which 
generally  is  the  average  temperature.  This  average  temperature  can  be  compared  with  the 
maximum  and  minimum  temperatures  of  the  components  from  the  FEA  simulation.  The 
numbers  shown  in  the  Figure  7  are  the  same  as  those  in  the  lumped  node  network. 

Table  5-  The  Analogy  of  the  Equivalent  Thermal  Circuit 


Electrical  Circuit 

Thermal  Circuit 

V 

[V] 

Voltage 

T 

[Cl 

Temperature 

I 

[AJ 

Current  Flow 

Q 

[WJ 

Heat  Flow 

a 

[1/an] 

Electrical 

Conductivity 

k 

[wrcm] 

Thermal 

Conductivity 

R 

[Q] 

Resistance 

Re 

[°C/W] 

Thermal 

Resistance 

C 

[F] 

Capacitance 

Ce 

[J/°C] 

Thermal 

Capacitance 

conductive  — 

k*Ac 

1 

convective  — 

h*  As 

Ce  =p*V*Cp 


(26) 


(27) 


where 

L  and  A  c  are  length  and  cross-section  area  of  the  components 
V  and  As  are  volume  and  surface  area  of  the  components 
p  and  CP  are  density  and  specific  heat  of  the  components 

Eq.  26  and  Eq.  27  show  how  to  calculate  the  R  and  C  values  for  a  simple  one-dimensional 
geometry.  Accurate  estimates  for  R  and  C  values  are  not  possible  for  complicated  geometries 
such  as  an  electric  motor.  A  detailed  3-D  solid  model  of  a  target  motor  is  constructed  and 
imported  to  FEA  software  like  ANSYS  to  perform  a  steady-state  and  transient  simulation.  With 
these  results  we  can  estimate  and  select  the  values  for  the  thermal  resistances  and  capacitances. 
The  advantage  of  this  method  is  that  we  can  evaluate  the  fidelity  of  the  lumped  node  model  with 
a  real  motor,  add  or  reduce  nodes  to  increase  the  model  accuracy  and  efficiency,  and  estimate  the 
maximum  error  between  the  node  temperature  and  maximum  temperature  in  a  real  motor. 

This  procedure  can  also  be  extended  to  model  the  gear-box,  motor-driver  and  drive-train,  and 
even  include  the  aircraft  wing  surfaces  and  frames.  The  thermal  network  can  also  be 
incorporated  into  a  multi-physics  model  to  simulate  the  electrical,  thermal,  and  mechanical 
performance  of  the  whole  EMA  and  its  supporting  structure.  Once  the  proper  thermal 
resistances  and  capacitances  are  selected,  this  simulation  engine  can  be  used  with  various  time 
dependent  boundary  conditions  during  the  whole  mission  duration  including  the  air  temperature, 
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aircraft  speed,  altitude  and  sunlight.  The  computational  requirement  of  such  a  simulation  is 
negligible  compared  to  the  FEA  simulation. 

The  node  network  is  constructed  as  shown  in  Figure  11.  The  electric  motor  has  two  thermal 
paths  with  other  parts  of  EMA  and  environment.  The  first  path,  also  the  major  one,  is  through 
the  convection  resistances  R5,  R5a  and  R5b,  and  the  radiation  resistance  R5r  (which  will  be 
discussed  later  in  this  chapter).  The  convection  resistances  represent  the  thermal  resistances 
between  motor  surface  and  the  air  around  it.  Considering  that  the  motor  will  be  installed  in  a 
small  space  inside  the  wing  structure,  this  convection  is  likely  to  be  natural  convection  or  weak 
forced  convection.  In  ANSYS  we  assume  the  heat  transfer  coefficient  for  the  forced  convection 
to  be  h=100  W/m2K,  which  is  a  typical  value  for  forced  convection  in  air  (Figure  9).  The  value 
of  R5,  R5a  and  R5b  can  be  easily  changed  to  accommodate  the  actual  flow  situation  in  the 
motor’s  final  installation  and  flow  conditions.  In  this  model,  the  epoxy  around  the  copper 
winding  is  combined  with  the  copper  windings  to  form  node  1 ,  due  to  the  very  small  temperature 
difference  between  them.  In  some  types  of  motors,  such  as  those  with  long  end  turns,  the  epoxy 
temperature  can  be  quite  different  than  that  of  the  copper  winding.  Separating  the  epoxy  and 
copper  winding  provides  more  flexibility  to  the  lumped-node  model,  and  allows  the  model  to  be 
applicable  to  different  motor  designs. 


Detailed  Lumped  node  model  Of  EMA  (transient) 


1.  copper 

la.  copper  in  front  end  turn 

lb.  copper  in  back  end  turn 

2.  epoxy  in  Iron 

2a.  front  end  turn  epoxy 
2b.  back  end  turn  epoxy 

3.  stator 

4.  case 

5.  case  surface 

5a.  front  cover  surface 
5b.  back  cover  surface 

6.  shaft  body 

7.  shaft  front  end 

8.  shaft  back  end 

9.  magnet 

10.  air  gap 

11.  connection  gear 

12.  front  bearing 

13.  back  bearing 

14.  front  cover 

15.  back  cover 

16.  wing  surface 

17.  gear  box  case 

18.  moving  mechanism 


R18  Surface  R16 


Gear  mechanism  assembly 


Figure  11  -  Lumped  Node  Model  of  EMA  (Transient) 


The  second  thermal  path  is  the  conductive  heat  transfer  from  the  motor  rotor  to  the  gear 
mechanism.  This  path,  via  R7a,  is  a  minor  heat  flow  path  because  the  rotor  and  bearing  loss  is 
relatively  small.  But  there  is  a  possibility  that  the  gear  mechanism  has  higher  temperature  than 
the  rotor  and  conduct  heat  back  to  motor.  There  could  be  additional  thermal  paths,  such  as  from 
the  motor  back  plate  to  the  aircraft  frame,  which  are  not  shown  in  Figure  11.  However, 
incorporating  additional  thermal  paths  into  the  network  and  changing  the  lumped-node  model 
accordingly  is  straightforward  with  the  method  presented  in  this  chapter. 
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The  following  procedure  is  used  to  determine  the  R  and  C  values  in  the  lumped  node  network  of 
motor  assembly  shown  in  Figure  1 1 .  The  thermal  resistances  and  capacitances  are  initially 
estimated  using  by  Eq.  26  and  Eq.  27.  The  R  values  can  be  adjusted  and  updated  by  comparing 
with  lumped  node  model  results  with  the  steady-state  results  obtained  from  ANSYS.  Table  4 
lists  the  nodal  temperatures  from  ANSYS  steady  state  simulation  for  motor  with  forced 
convection  boundary  conditions.  After  the  R  values  have  been  determined,  the  thermal 
capacitances  can  be  updated  by  comparing  the  transient  period  of  the  lumped  node  results  with 
the  transient  temperature  distribution  obtained  from  ANSYS.  Table  6  lists  the  R  and  C  values 
used  in  the  lumped  node  model. 

Table  6  -  The  R  and  C  Values  of  Lumped  Node  Model  of  the  Motor 


R 

(K/W) 

C 

(J/K) 

R1 

0.0029 

Cl 

77.16 

Rla 

0.115 

Cla 

25.94 

Rib 

0.115 

Clb 

25.94 

Rlc 

0.896 

C2 

11.95 

Rid 

0.896 

C2a 

44.69 

R2 

0.0093 

C2b 

44.69 

R2a 

0.395 

C3 

291.21 

R2b 

0.395 

C4 

499.39 

R2c 

1.045 

C6 

232 

R2d 

1.045 

C7 

35.58 

R3 

0.0305 

C8 

24.51 

R5  (h=  1 00) 

0.492 

C9 

19.43 

R5  (h=  1 0) 

4.92 

C12 

140.17 

R5a  (h=100) 

3.199 

C13 

11.23 

R5a  (h=10) 

31.99 

C14 

186.65 

R5b  (h=100) 

3.199 

C15 

186.65 

R5b  (h=10) 

31.99 

R6 

0.03556 

R6a 

4.054 

R6b 

3.22 

R7 

1.961 

R8 

1.592 

R9 

9.21 

RIO 

19.76 

R13 

0.86229 

R14 

0.76784 

R14a 

0.254 

R15a 

0.254 

R15b 

0.18 
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After  all  the  resistor  and  capacitor  values  of  the  lumped-node  network  model  of  the  motor  in 
Figure  1 1  are  known,  the  same  values  of  the  resistors  and  capacitors  can  be  used  to  simulate  the 
temperature  response  of  the  motor  parts  with  any  combination  of  heat  losses  and  boundary 
conditions.  For  this  standard  resistor-capacitor  electrical  network,  a  set  of  first-order  ordinary 
differential  equations  can  be  written  as  (Eq.  28  -  Eq.  44): 
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where 

Tamb  is  ambient  temperature, 

T bay  is  bay  wall  temperature, 

0  is  heat  generation 

This  equation  set  can  be  solved  by  standard  numerical  methods.  In  Table  4,  the  temperatures 
with  the  presence  of  radiation  heat  transfer  is  also  considered  and  listed  in  the  last  two  columns. 

It  should  be  noted  that  in  a  situation  where  air  cooling  is  not  effective  (such  as  natural  convection 
with  a  low  h),  the  motor  temperature  could  be  very  high  and  radiative  heat  transfer  becomes  the 
dominant  mode  of  heat  transfer.  As  radiation  could  be  important,  this  feature  is  included  in  Eq. 
10.  The  parallel  radiation  thermal  path  is  modeled  by  adding  a  radiative  thermal  resistance  to  the 
lumped  node  model.  The  radiative  thermal  resistance  value  can  be  obtained  by  using  the 
following  equations. 


K  =  Eo(r12  +  r22)(T1  +  7'2) 


(45) 


77.  = 


h.A 


(46) 


where 

e  is  the  emissivity 

cris  the  Stefan-Boltzmann  constant 

T i  is  motor  casing  temperature 

T 2  is  bay  wall  temperature 

A  is  surface  area 

Since  the  motor  casing  temperature  is  a  temperature  to  be  determined,  the  inclusion  of  radiation 
contribution  involves  an  iterative  process  in  the  solution  procedure  to  solve  Eq.  28  -  Eq.  44. 

3.8.2. 1  Examples 

The  simulation  results  of  the  lumped  node  model  are  compared  with  the  FEA  results  under  the 
same  initial  condition  and  boundary  conditions  (the  step  heat  load  given  in  Figure  8,  Q  =  600  W, 
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2 

h=100  W/m  K).  Figure  12  shows  the  lumped-node  model  and  the  FEA  produce  essentially  the 
same  temperature  for  the  copper  winding. 


Figure  12  -  Lumped-Node  Network  Simulation  Result 

This  lumped  node  model  is  also  tested  with  a  pulsed  heat  load,  where  the  thermal  load  is  a  square 
wave  form  repeating  with  time,  as  shown  in  Figure  13.  Figure  14  shows  a  comparison  between 
the  ANSYS  results  with  the  lumped  node  model  for  the  copper  winding  temperature.  The 
excellent  agreement  between  ANSYS  results  and  lumped-node  model  proves  that  this  model  is  a 
very  effective  and  attractive  alternative  to  the  computational  intensive  FEA  simulation.  It  is 
much  faster  to  obtain  essentially  the  same  results  as  the  FEA  by  using  the  lumped  node  model 
(less  than  a  minute  for  the  lumped  node  model  compared  to  more  than  20  hours  for  the  FEA  on 
same  computer  to  simulate  750  seconds  of  heat  load  profile). 


Figure  13  -  Pulsed  Heat  Load 
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Figure  14  -  Winding  Temperature  Comparison  (FEA  and  Lumped  Node) 

A  second  time  varying  heat  load  is  simulated  to  show  the  capability  of  the  lumped  node  model. 
Shown  in  Figure  15  is  a  highly  dynamic  heat  load  typical  of  those  of  a  10  hp  motor  used  for 
powering  a  flight  control  surface.  The  transient  temperature  of  the  copper  winding  and  the  motor 
case  are  given  in  Figure  16. 
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Figure  15  -  Heat  Load 
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Figure  16  -  Winding  and  Motor  Case  Temperatures 

3.8.3  Lumped  Node  Model  with  PCM  Module 

In  situations  where  the  heat  load  is  high  and  the  convection  heat  transfer  is  ineffective  (such  as  in 
an  aircraft  bay  with  near  stagnant  air),  the  temperature  of  the  motor  components  can  reach  an 
unacceptable  level.  For  example,  for  a  steady  heat  load  of  600  W  with  h=10  W/m  K,  the  copper 
winding  temperature  reaches  nearly  400°C  (see  last  column  of  Table  4).  Another  situation  where 
excessively  high  temperature  can  occur  in  the  copper  winding  even  with  forced  convection 
(h=100  W/m  K)  is  when  the  ambient  temperature  is  high.  With  an  ambient  temperature  of 
100°C  instead  of  22°C,  the  copper  winding  with  h=100  W/m2K  could  reach  200°C.  A  solution 
to  keep  the  motor  temperature  at  an  acceptable  level  is  to  use  phase  change  materials  (PCM)  with 
a  melting  temperature  below  150°C  to  store  the  excess  thermal  energy.  Figure  17  shows  the  3-D 
model  of  the  motor  and  an  aluminum/PCM  composite. 


0.030 


Figure  17  -  Aluminum/PCM  Composite 
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We  can  embed  phase  change  materials  in  an  aluminum  foam  structure  to  form  an 
aluminum/phase  change  material  (PCM)  composite.  To  keep  the  copper  winding  below  150°C 
by  storing  the  excess  thermal  energy  in  the  PCM,  we  should  choose  a  PCM  with  a  melting 
temperature  range  below  150°C  and  a  high  latent  heat  of  fusion.  In  this  chapter,  we  select 
Erythritol  as  the  PCM  [21].  This  PCM  has  a  high  latent  heat  (340  kJ/kg)  and  suitable  phase 
change  temperature  range  (~120°C).  The  solid  and  liquid  specific  heats  are  1.18  kJ/kg  and  2. 1 5 
kJ/kg,  respectively.  The  PCM  is  modeled  as  a  material  with  a  temperature-dependent  specific 
heat  as  shown  in  Figure  18. 


Figure  18  -  Specific  Heat  vs.  Temperature 

The  lumped-node  network  shown  in  Figure  19  is  based  on  the  motor  assembly  geometry. 
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Figure  19  -  Lumped-Node  Model  Including  PCM 
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The  additional  R  values  in  this  network  are  listed  in  Table  7  where  R5,  R5a  and  R5b  are  thermal 
contact  resistances  between  the  case,  front  cover  and  back  cover  to  the  PCM  module.  The  PCM 
thermal  capacitance  is  temperature  dependent  and  is  obtained  by  multiplying  the  specific  heat 
from  Figure  18  with  the  mass  of  the  PCM  (see  Eq.  47). 

Table  7  -  Network  Parameters  for  Equivalent  Thermal  Resistances 


R 

(K/W) 

R5 

0.0018 

R5a 

0.01 

R5b 

0.01 

R19  (h=10) 

2.72 

R19  (h=100) 

0.34 

To  solve  the  PCM  behavior,  Eq.  47  is  added  into  the  lumped  node  model,  and  Eq.  35,  43  and  44 
are  replaced  by  Eq.  35’,  43 ’and  44’. 
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3.8.3. 1  Examples 

Two  heat  load  cases  are  simulated  using  the  lumped  node  thermal  network  model.  The  cases  are 
listed  in  Table  8.  In  both  of  these  cases,  the  ambient  temperature  was  chosen  to  be  100°C. 

Table  8  -  Lumped-Node  Cases 


Case 

No. 

h 

(W/m2-K) 

Load 

Type 

Load  Power( 
W) 

1 

10 

step 

600 

2 

100 

step 

600 

Figures  20  and  21  show  the  simulation  results  for  the  copper  winding  temperature.  The  results 
show  that  PCM  module  can  hold  the  copper  winding  temperature  under  130°C  for  as  long  as  400 
s  (h=10  W/m2-K)  or  500  s  (h=100  W/m2-K)  under  the  600  W  heat  load. 
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Figure  20  -  h=10  W/m2K,  PCM  (Case  1) 


Figure  21  -  h=100  W/m2K,  PCM  (Case  2) 

The  lumped  node  thermal  network  is  a  feasible  approach  to  simulate  the  mission-level  EMA 
thermal  performance.  A  new  lumped  node  modeling  technique  which  is  based  on  tuning  the  R 
and  C  values  with  the  FEA  results  is  developed  to  overcome  the  difficulty  in  obtaining  accurate 
R  and  C  values  in  the  motor  thermal  network.  The  lumped  node  model  has  been  successfully 
applied  to  step  and  transient  heat  loads  with  cooling  by  both  convection  and  radiation.  For 
situations  where  phase  change  materials  may  be  needed  to  store  the  thermal  energy  to  avoid  the 
motor  reaching  excessively  high  temperature,  a  PCM  module  has  been  incorporated  successfully 
into  the  lumped  node  model  as  well. 

3.9  An  EMA  Example  and  Results 

A  commercially  available  EMA  was  measured,  modeled,  simulated,  and  experimentally  tested. 
Here  are  the  results. 
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a.  Parameters 


The  key  parameters  for  the  test  EM  A  are  shown  in  Table  9. 

Table  9  -  EMA  Parameters 


Symbol 

Name 

Value 

Ncr 

Gearing  ratio 

1963  rad/m 

s 

Number  of  slots 

12 

p 

Number  of  poles 

10 

Ff 

Power  train  friction 

342  N 

m 

Rod  mass 

8.5  kg 

I 

Rotor  moment  of  inertia 

1 13.2  kg-mm2/rad2 

R s,ref 

Phase  resistance  at  t  , 

1.4  Q 

Tref 

Reference  temperature 

20  °C 

CCr 

Temperature  coefficient  of  rs 

0.004041  1/°C  [22] 

^ PM 

PM  flux  linkage 

0.149  Wb 

^ d,peak 

Peak  direct  inductance 

17.35  mH 

^ q,peak 

Peak  quadrature  inductance 

17.27  mH 

Figure  22  shows  the  force-stroke  data  for  the  friction  analysis. 
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Figure  22  -  Static  Friction  Analysis  Showing  Stroke  (top)  and  Load  Force  (Bottom).  The 
Green  Lines  on  the  Force  Plot  Represent  the  Mean  Forces  to  Pull  (-258  N)  or  Push  (425  N) 
the  EMA  Rod,  Overcoming  the  Friction  of  the  Drive  Train  (342  N).  The  Bias  towards 
Force  Required  to  Push  is  the  Weight  of  the  EMA  Rod  (84  N) 

The  direct  and  quadrature  inductances  are  obtained  as  shown  in  Figure  23  and  Figure  24.  The 
FEM  model  is  simulated  with  a  range  of  currents  0  A  to  50  A,  a  range  of  torque  angles  0  to  n 
radians,  and  a  range  of  rotor  angles  0  to  n/3  radians.  For  this  motor,  a  range  of  n/3  radians 
captures  one  full  cycle  of  variation  in  inductance  values  due  to  rotor  angle. 


Figure  23  -  Direct- Axis  Inductance  as  a  Function  of  Direct  Current, /,, ,  and  Quadrature 

Current,  iq 
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Figure  24  -  Quadrature-Axis  Inductance  as  a  Function  of  Direct  Current,^,  and 

Quadrature  Current,  iq 

A  computer  rendering  of  the  motor  geometry  is  shown  in  Figure  25.  This  geometry  is  used  for 
FEM  modeling. 


Figure  25  -  Motor  Geometry  (Left)  and  Slot  Measurements  (Right) 

Table  10  shows  some  of  the  material  properties  of  the  EMA  motor. 

Table  10  -  EMA  Material  Properties 


Part 

Name 

cP 

k 

P 

Casing 

AL-384 

963 

96.2 

2820 

Lamination 

Steel  1010 

448 

51.9 

7872 

Windings 

Copper 

385 

401 

8940 

Epoxy 

(unknown) 

1050 

0.85 

1200 

Axle 

Steel 

480 

15 

7750 

Magnets 

NdFeB-30 

5e4 

8.95 

7400 

Cp  is  the  specific  heat  capacity. 
k  is  the  thermal  conductivity. 
p  is  the  density. 
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b.  Results 


The  desired  stroke  profile  and  the  motor's  actual  stroke  profile  are  shown  in  Figure  26.  The 
motor  followed  the  desired  stroke  profile  very  well. 


Stroke  Comparison 


Figure  26  -  Desired  (Blue)  and  Actual  (Green)  Stroke  Profiles 


The  load  force  profile  and  the  motor-generated  force,  are  shown  in  Figure  27. 


-15o 


3 

Time,  t  [s] 


Figure  27  -  Input  Load  Force  Profile  (Blue)  and  Motor-Generated  Force  (Green).  The 

Load  Force  was  a  Sustained  450  N 

The  DQO  voltages  and  currents  are  shown  in  Figure  28  and  Figure  29,  respectively. 
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Figure  28  -  DQO  Voltages  with  Pulse  Form  due  to  PWM 


Currents 


Figure  29  -  DQO  Currents 

If  the  motor  is  treated  as  a  system,  the  powers  into  (positive)  and  out  of  (negative)  this  system 
are  shown  in  Figure  30. 
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Terminal  power 


Time,  t  [s] 

Motor-generated  power 


Figure  30  -  Motor  Powers  Showing  the  Total  Electrical  Power  at  the  Terminals  of  the 
Motor  (First),  the  Power  Loss  in  the  Windings  (Second),  the  Power  Stored  in  the  Magnetic 
Fields  (Third),  and  the  Mechanical  Power  Output  on  the  Rotor  (last) 

The  electrical  power  at  the  terminals  is  the  power  on  the  DC  bus  line.  Note  that  sometimes  this 
power  is  negative,  meaning  that  the  motor  is  actually  sending  power  back  to  the  source.  The 
reason  for  this  is  that  the  wind  load  is  actually  driving  the  flight  control  surface,  pushing  against 
it  as  it  yields  to  the  wind  pressure. 

If  the  drive  train  is  treated  as  a  system,  the  power  into  and  out  of  this  system  is  shown  in  Figure 
31. 
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Motor-generated  power 


Time,  t  [s] 


Figure  31  -  Drive  Train  Powers  Showing  the  Mechanical  Power  Input  on  the  Rotor 
(First),  the  Load  Power  (Second),  the  Frictional  Power  (Third),  the  Gravitational  Power 

(Forth),  and  the  Net  Kinetic  Power  (Last) 
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The  comparison  of  the  recorded  temperatures  and  the  simulated  temperature  of  the  stator  is 
shown  in  Figure  32. 


Temperature  Comparison 


Figure  32  -  Temperature  Comparison  (Above)  and  Again  (Below)  Zoomed  in.  These  Plots 
have  Environmental  Temperature  (blue),  experimentally  Measured  Temperature  of  Stator 
(Green),  Simulated  Temperature  of  Stator  (Red),  and  the  Simulated  Temperature  of  the 

Windings  (Cyan) 

The  difference  here  is  actually  quite  small,  only  about  0.34  °C,  much  less  than  the  1  °C 
sensitivity  of  the  thermocouples  used. 

The  accuracy  of  the  temperature  prediction  is  more  readily  observed  for  longer  simulation 
profiles.  Another  test  is  shown  in  Figure  33  lasting  a  half-hour.  The  simulation  of  this  test 
matches  the  stator  temperature  very  well  with  a  maximum  error  of  0.17  °C. 
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Temperature  Comparison 


Figure  33  -  Temperature  Comparison  for  a  Half-Hour  Run  Showing  Environmental 
Temperature  (Blue),  Experimentally  Measured  Temperature  of  Stator  (Green),  Simulated 
Temperature  of  Stator  (Red),  and  the  Simulated  Temperature  of  the  Windings  (Cyan). 
Note  that  the  Simulated  Temperatures  of  the  Windings  and  Stator  are  About  the  Same, 
Indicating  the  Low  Thermal  Resistance  in  the  Motor 

The  test  shown  in  Figure  32  is  run  using  kv  and  ka  control  coefficient  values  of  0.4  and  0.85, 

respectively.  When  these  coefficients  are  instead  both  set  to  1.04,  the  mean  copper  loss  in  the 
windings  increased  by  a  factor  of  1 1.5.  The  temperature  comparison  between  simulation  and  test 
is  shown  in  Figure  34  and  the  terminal  power  is  shown  in  Figure  35. 


Temperature  Comparison 


Figure  34  -  Temperature  Comparison  Showing  the  Disparity  between  the  Experimentally 
Measured  Temperature  of  the  Stator  (green)  and  the  Simulated  Temperature  (red)  when 
Different  Values  of  Control  Coefficients  are  used.  In  this  Case,  the  Test  Data  is  shown  Only 
for  Reference.  Obviously,  with  this  Control  Tuning,  the  Temperatures  are  not  expected  to 

Match 
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Figure  35  -  Peak  Terminal  Power.  Note  that  the  Peak  Power  here  is  Over  50%  Higher 

than  Previously  in  Figure  31 

This  illustrates  the  significance  that  control  has  on  excess  heat  generation.  Inadequate  control  can 
lead  to  significantly  higher  losses. 

There  is  a  strong  motivation  to  move  towards  all-electric  aircraft.  However,  the  significant 
obstacles  to  full  realization  of  all-electric  aircraft  are  the  thermal  management  of  heat  loads  and 
power  management.  In  this  chapter,  it  shows  the  successful  implementation  of  an  EMA  model 
with  a  control  algorithm,  a  PWM  component,  an  electromechanical  component  with  nonlinear 
parameters,  and  a  thermal  component.  The  model  takes  a  desired  position  profile,  a  load  force 
profile,  and  an  environmental  temperature  profile  as  inputs  and  predicts  the  thermal  behavior  of 
the  motor  and  the  transient  power  demands  of  the  EMA.  It  is  shown  that  control  plays  a 
powerful  role  in  how  much  excess  heat  is  generated  and  how  much  the  power  demand  spikes.  It 
is  also  shown  that  the  temperatures  within  the  motor  of  the  EMA  are  all  within  a  fairly  narrow 
band.  This  approach  of  using  finite  element  method  and  nonlinear,  lumped-element,  integrated 
modeling  gives  a  complete  picture  of  the  dynamic  behavior  of  an  EMA,  thereby  helping  to  make 
all-electric  aircraft  a  reality. 
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4.  MODELING  AND  SIMULATION  OF  POWER  LOSS  IN 
DRIVER  UNIT  OF  ELECTROMECHANICAL  ACTUATOR 

4.1  Introduction 

Since  the  operation  of  the  EMA  is  highly  transient  and  unpredictable  due  to  the  unexpected 
turbulence,  the  load  torque  and  stroke  might  change  rapidly  or  need  to  keep  stable  for  a  certain 
amount  of  time.  The  power  loss  and  heat  generation  is  different  and  more  challenging  from 
ordinary  motors  that  operates  under  normal  conditions  with  torque  and  speed  that  does  not 
change  dramatically.  The  highly  transient  load  and  stroke  profile  will  result  in  a  highly  transient 
current  and  voltage  in  the  driver  circuit.  And  there  may  be  a  large  amount  of  regenerative  power 
generated  from  the  deceleration  process  of  the  EMAs,  which  may  also  be  converted  to  heat. 

In  this  chapter,  we  study  the  heat  problem  of  the  control  and  drive  units  of  EMAs,  and  build  a 
model  to  calculate  and  simulate  the  power  loss  and  heat  generation  in  the  driver  board.  The 
driver  unit  consists  of  a  power  inverter,  a  power  dissipating  resistor  and  a  control  circuit.  The 
power  loss  of  each  part  is  studied.  The  heat  loss  in  the  power  inverter  comes  mainly  from  the 
power  switches  -  IGBTs.  The  on-state  loss  is  proportional  to  the  current  of  the  motor,  and  the 
switching  loss  is  determined  by  the  switching  frequency  as  well  as  currents  and  voltages.  The 
power  loss  in  the  power  dissipating  resistor  is  determined  by  the  regenerative  power,  the 
capacitor  and  the  control  algorithm  to  stabilize  the  bus  voltage,  which  varies  from  different 
mission  profiles  and  different  applications.  All  those  parameters  can  be  obtained  in  our 
simulation  code.  The  power  loss  on  the  control  circuit  is  negligible  compared  with  the  power  loss 
on  IGBTs  and  the  power  dissipating  resistor,  and  generates  very  little  heat  in  the  system. 

4.2  PWM  Drive  Simulation 
a.  Background 

Pulse-width  modulation  (PWM)  is  a  commonly  used  technique  for  controlling  power  in  inertial 
electrical  devices,  made  practical  by  modem  electronic  power  switches. 

The  average  value  of  voltage  (and  current)  fed  to  the  load  is  controlled  by  turning  the  switch 
between  supply  and  load  on  and  off  quickly.  The  longer  the  switch  is  on  compared  to  the  off 
periods,  the  higher  the  power  supplied  to  the  load. 

The  main  advantage  of  PWM  is  its  very  low  power  loss  in  the  switching  devices.  When  a  switch 
is  off,  there  is  practically  no  current,  and  when  it  is  on,  there  is  almost  no  voltage  drop  across  the 
switch.  Power  loss,  being  the  product  of  voltage  and  current,  is  thus  in  both  cases  close  to  zero. 
PWM  also  works  well  with  digital  controls,  which,  because  of  their  on/off  nature,  can  easily  set 
the  needed  duty  cycle. 
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Carrier  Wave 


The  PWM  wave  generation  is  shown  in  Figure  36.  Comparing  the  carrier  wave  and  modulating 
wave  determines  the  PWM  output.  In  this  example,  there  are  two  output  levels  of  the  PWM 
wave,  high  and  low.  The  amplitude  of  the  carrier  or  sinusoidal  wave  directly  corresponds  to  the 
output  level:  If  the  carrier  wave  is  higher  than  the  modulating  wave,  the  output  level  is  high; 
otherwise,  the  output  level  is  low.  Correct  configuration  of  the  modulating  wave  and  output 
voltage  can  create  the  output  with  two  voltages  which  can  be  used  to  drive  the  motor.  The  same 
result  occurs  because  the  main  component  in  the  square  wave  is  still  coordinated  with  the  carrier 
wave.  The  high-  frequency  components  that  made  it  square  would  be  automatically  filtered  by 
the  inductance  and  would  not  contribute  to  the  power  input  to  the  motor.  This  attribute  allows 
us  to  use  the  PWM  wave  to  drive  the  motor. 

b.  Simulation  Model 

The  model  used  for  simulation  consists  of  a  PMSM  motor  and  a  controller  previously  developed 
for  this  project.  Two  inputs  are  given:  a  time-history  profile  (a  sequence  of  data  points  versus 
time)  of  desired  positions  (stroke  reference  profile)  and  another  time-history  profile  of  load 
forces  (force  profile).  Both  of  these  profiles  can  be  translated  to  the  perspective  of  the  motor 
through  the  conversion  factor  of  the  drive  train  ratio.  This  enables  us  to  generate  a  rotor-angle 
profile  and  a  load-torque  profile.  Then,  for  each  time  step,  the  appropriate  applied  voltages  to 
make  the  motor  follow  the  stroke  profile  in  the  presence  of  the  load  torque  profile  are  calculated 
by  the  control  algorithm.  These  voltages  are  then  applied  to  the  motor,  at  which  point  the 
simulation  enters  a  convergence  loop  to  determine  the  correct  currents  and  motor-generated 
torques  in  the  presence  of  such  applied  voltages.  With  the  motor-generated  torque  determined, 
the  new  rotor  position  can  be  calculated  for  that  time  step.  Thus,  the  motion  of  the  motor  is 
simulated.  Now,  for  each  time  step,  once  the  currents  are  calculated,  the  power  loss  in  the 
windings  can  also  be  found.  These  copper  power  losses  then  become  the  heat  sources  for  the 
thermal  algorithm,  where  temperatures  for  various  points  (nodes)  in  the  motor  are  calculated. 
The  outputs  from  the  simulation  are  the  actual  stroke  profile  of  the  motor  (as  opposed  to  what 
was  desired:  the  reference  profile),  the  electrical  power  flowing  in  and  out  (regenerative)  of  the 
motor,  and  the  temperature  profiles  for  the  various  nodes. 

The  Dynamic  equations  of  the  motor: 
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(1  /  Rd  X  RJd  +  ^meLqiq  +  UdPWM  ) 
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dt 
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dt 
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dt 

where  ud  is  input  direct  voltage,  uq  is  input  quadrature  voltage,  id  is  direct  current,  i  is 
quadrature  current,  Rs  is  phase  resistance,  Ld  is  direct  inductance,  Lq  is  quadrature  inductance, 
(Ome  is  the  mechanical  frequency  multiplied  by  the  number  of  pole  pairs  p/2,  XPM  is  the  flux 
linkage  from  the  permanent  magnet,  J  is  the  rotor’s  moment  of  inertia,  tm  is  the  EM  torque,  and 
tl  is  the  load  torque,  0m  is  the  rotor’s  position,  com  is  the  mechanical  frequency. 


The  controller  follows  the  voltage  control  equation: 

(52) 

(53) 

(54) 

where  the  variables  with  stars  are  the  calculated  values  of  next  time  step. 

Besides  the  controller  and  the  motor,  a  PWM  block  is  introduced.  The  input  to  this  PWM  block 
is  the  control  voltage,  ud  and  uq .  The  output  is  the  PWM  control  voltage.  This  module  adjusts  the 

control  signal,  generates  a  PWM  wave,  and  sends  it  to  the  motor  model.  The  whole  system’s 
structure  is  shown  in  Figure  37. 
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Figure  37  -  Simulation  System  Structure 


c.  PWM  Drive 


The  PWM  drive  module  is  added  to  our  simulation  to  realize  PWM  function.  The  input  is  the 
control  signal,  ud  and  aq . 

Before  modulation,  the  input  signal  needs  to  be  converted  by  the  Park  Transformation.  This  is  a 
mathematical  conversion  that  changes  the  quantities  between  stator-referenced  domain  (abc 
frame)  into  rotor-referenced  domain  (dqO  frame).  In  a  simulation,  we  use  a  dqO  frame  because  it 
simplifies  the  computation.  However,  modulation  of  the  voltages  should  be  done  in  an  abc  frame 
because  that  is  the  quantities  with  physical  significance  and  can  be  detected  and  used  in  real 
applications. 

The  basic  equation  of  Park  Transformation  is: 
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where  0m  is  rotor  mechanical  angle 

In  the  PWM  drive  module,  the  quantities  to  be  transformed  are  ud  and  aq .  So  the  Park 
Transformation  can  be  written  as: 
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PWM  is  done  by  comparing  the  carrier  wave  with  the  control  signal.  The  carrier  waves  are  often 
high-frequency  period  signals.  The  carrier  wave  frequency  carrier  wave  should  be  higher  than 
the  highest  frequency  component  that  we  want  to  keep  in  the  control  signal.  That  is  because, 
after  modulation,  the  components  with  frequency  higher  than  the  carrier  wave  will  be  lost. 


The  time  step  is  an  important  parameter  in  the  simulation.  After  PWM  is  introduced  into  our 
simulation,  the  time  step  should  be  smaller  because  of  the  high-frequency  carrier  wave’s 
introduction.  This  will  significantly  increase  the  simulation  time.  Some  strategy  is  used  to  reduce 
the  computation  time.  The  simulation  uses  small  time  steps  when  doing  computation  with  PWM 
but  for  other  parts  a  lower  resolution  is  used. 


The  PWM-driven  simulation  results  show  that  our  motor  and  controller  model  can  function 
properly  with  PWM  driving  the  motor.  Some  simulation  results  and  comparison  with  simulation 
results  without  PWM  are  shown  in  Figure  38  and  Figure  39: 
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Figure  38  -  Mechanical  Stroke  Following  Results  Comparison 
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Figure  39  -  D-Axis  and  Q-Axis  Voltages  Comparison 
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4.3  Structure  of  the  Electronic  Unit 

The  electronic  unit  is  a  circuit  board  combining  the  control  and  drive  modules  for  the  EMA.  It  is 
based  on  a  Danaher  control  unit  (SERVOSTAR  S610).  The  control  module  has  a  microprocessor 
and  its  peripheral  circuits.  The  drive  module  is  an  IGBT  package.  There  is  a  resistor  called  an 
unloading  resistor.  Its  function  is  to  dissipate  the  excessive  regenerative  energy  produced.  There 
is  a  capacitor  on  the  DC  bus  to  store  the  regenerative  power.  The  amount  of  energy  that  the 
capacitor  could  store  is  limited  by  the  DC  bus  voltage  range.  When  the  DC  bus  voltage  reaches 
its  upper  limit,  the  capacitor  cannot  be  charged,  and  the  excessive  regenerative  energy  will  be 
dissipated  on  the  resistor.  The  circuit  board  is  shown  in  Figure  2.  The  red  square  indicates  where 
the  IGBTs  and  the  resistor  are.  The  IGBT  package  is  under  the  heat  spreader  marked  by  “EP”. 
The  green  resistor  is  the  one  for  dissipating  regenerative  power.  On  the  left  is  the  heat  sink  which 
covers  the  red  square  area.  It  is  removed  to  reveal  the  IGBT  package  and  the  unloading  resistor. 


Figure  40  -  Electronic  Unit  of  the  EMA 

The  board  can  be  divided  into  three  main  components:  controller  circuit,  the  IGBT  package  and 
unloading  resistor.  The  diagram  that  shows  the  function  and  relationship  of  each  component  is  in 
Figure  41.  The  bold  wires  indicate  the  power  flow  through  and  the  current  can  be  large.  The  thin 
wires  means  those  electrical  quantities  are  just  signals  and  do  not  carry  large  current. 
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Figure  41  -  Electronic  Unit  Block  Diagram 

4.4  Power  Loss  on  the  EU 

Since  we  are  concerned  about  power  loss  that  will  cause  heating  problems  and  do  not  focus  on 
efficiency,  only  the  absolute  amount  of  power  loss  is  meaningful  to  our  study. 

As  shown  in  Figure  41,  large  power  only  goes  through  the  IGBTs,  the  unloading  resistor  and  the 
capacitor.  The  controller  module  only  deals  with  signals  and  does  not  carry  large  currents.  The 
power  level  of  that  part  is  much  smaller  than  the  other  two  parts  on  the  board. 

4. 4. 1  Power  Loss  on  IGBTs 

IGBTs  are  switching  components  that  control  the  input  voltage  to  the  motor.  In  this  EM  A,  a  full 
bridge  IGBT  power  module  is  used  as  an  inverter  to  drive  the  motor.  There  are  six  IGBTs 
packaged  in  the  module.  It  connects  to  the  DC  source  which  provides  the  power,  and  provides 
PWM  voltage  to  the  motor.  The  switches  are  controlled  by  the  controller,  and  operate  at 
approximately  1  kHz. 

In  an  aircraft  power  system,  a  DC  power  supply  is  available  to  the  EMA  system.  In  order  to  drive 
the  motor,  we  need  to  use  the  inverter  to  generate  AC  power.  An  inverter  is  also  essential  to 
implement  PWM  control.  A  typical  full  bridge  inverter  diagram  is  shown  in  Figure  42.  Six 
IGBTs  are  used  to  build  the  inverter. 
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Figure  42  -  Two  Level  Voltage  Source  Inverter,  SI  ~  S6  are  Switches,  Which  are  IGBTs 

in  the  EMA  Application 

The  IGBTs  are  controlled  by  the  gate  signals,  and  function  as  switches.  The  controller 
determines  when  the  IGBTs  should  be  turned  ON  or  OFF.  The  control  of  the  motor  is 
accomplished  by  switching  IGBTs.  The  motor  stator  current  will  go  through  the  IGBTs  that  are 
ON,  and  the  motor  terminal  voltage  will  be  applied  on  the  IGBTs  that  are  OFF.  The  voltage  and 
current  can  be  large  if  the  power  consumed  by  the  motor  is  high. 

Since  the  IGBTs  deal  with  large  current  and  voltage,  the  power  loss  on  them  are  very  important 
in  the  system.  The  power  loss  on  IGBTs  comes  from  two  sources:  on-state  loss  and  switching 
loss. 

On-state  loss  is  caused  by  the  current  flowing  through  the  non-perfect  device’s  on-state 
resistance.  While  the  device  is  turned  “ON”  and  current  flows  through  the  device,  there  will  be  a 
voltage  drop  and  the  on-state  resistance  could  be  used  to  characterize  this  power  loss. 

Switching  loss  the  loss  happens  when  the  IGBT  is  turning  on  and  off.  The  transition  between  ON 
and  OFF  states  is  not  instant,  during  the  rise  and  fall  of  the  current  and  voltage  the  product  of 
them  will  generate  some  power  loss. 

For  a  given  IGBT,  the  on-state  loss  is  related  to  the  current  and  voltage  on  the  IGBTs,  which 
depends  on  the  control  scheme  and  the  load  of  the  motor.  The  switching  loss  is  determined  by 
the  switching  frequency. 


48 

Approved  for  public  release;  distribution  unlimited. 


a.  On-state  Loss  of  the  IGBTs 


The  IGBTs  function  as  switches  in  the  circuit.  When  the  switch  is  turned  “ON”,  it  should 
perform  like  an  ideal  switch,  and  have  zero  resistance.  However,  for  IGBTs,  especially  those 
with  high  power  ratings,  the  voltage  drop  is  not  negligible  even  when  they  are  turned  ON. 

Since  the  current  is  relatively  large  when  the  switch  is  ON,  with  a  voltage  drop,  the  power  loss 
PoN  =  VONl 


where  vON  is  the  on-state  voltage  drop  across  the  device,  and  i  is  the  current  flowing  through  the 
collector  (for  IGBTs)/drain  (for  MOSFETs)  terminal. 

Then  the  power  loss  can  be  expressed  as  follows,  with  RON  being  the  on-state  resistance: 


Pon  i  R 


ON 


(71) 


Figure  43  is  an  example  showing  the  current  versus  voltage  of  an  IGBT.  The  data  is  from  the 
datasheet  of  a  600Y/300A  IGBT  module. 


Collector  current  vs.  Collector-Emitter  voltage 
Tj=25°C 


Figure  43  -  Collector  Current  vs.  Collector-Emitter  voltage.  The  Curve  Becomes  Linear 
When  Current  is  higher  than  100 A.  In  that  Range,  an  Equivalent  Resistance  can  be  used  to 

calculate  the  Voltage  Drop 


As  shown  in  the  figure,  the  current  versus  voltage  is  almost  linear  when  current  is  higher  than 
100A.  An  equivalent  resistance  can  be  used  to  calculate  the  on-state  loss  if  we  know  the  stator 
current. 
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Because  the  stator  current  is  dependent  on  the  mission  profile,  the  on-state  loss  is  also 
determined  by  the  mission  profile.  Different  load  and  stroke  curves  may  result  in  very  different 
on-state  losses  on  the  IGBTs. 

b.  Switching  Loss  of  the  IGBTs 

When  the  device  is  turned  ON  or  OFF,  the  current  cannot  change  instantaneously,  there  is  a 
switching  time  for  it  to  rise/fall.  During  that  period  of  time,  both  current  and  voltage  are  not  zero 
or  negligible.  The  power  loss,  which  is  the  product  of  current  and  voltage,  needs  to  be 
considered.  The  switching  process  is  shown  in  Figure  44. 


WPUT  Signal 


Figure  44  -  Definition  of  Switching  Time  of  IGBT.  The  “Input  Signal”  is  the  Input  to  the 

Gate,  and  the  Device  is  an  N-Channel  IGBT 

Switching  loss  occurs  at  every  switching  of  the  device.  It  includes  on,  off,  and  reverse  recovery 
losses.  For  a  given  device,  the  power  loss  during  each  switch  operation  could  be  modeled  by  the 
equations  given  in  [23], 

E„  =  (72) 

where  Esw  is  the  switching  energy  loss  for  one  switching  operation,  i(t)  is  the  load  current, 

Asw  and  Bsw  are  curve  fitting  constants  for  the  device.  The  equation  can  be  used  for  turn-on, 
turn-off  losses  as  well  as  reverse  recovery  losses. 

Assuming  load  current  is  constant,  then  Esw  is  constant  for  a  specific  device,  and  the  switching 
loss  could  be  expressed  as: 

P,.=(E„  +  Evf  +  Err)-f  (73) 

where  /  is  the  switching  frequency.  Eon  is  turn-on  loss,  Eoff  is  turn-off  loss  and  Err  is  reverse 

recovery  loss.  It  shows  that  the  switching  loss  is  proportional  to  switching  frequency.  A  curve 
showing  switching  energy  loss  versus  current  is  usually  provided  in  the  datasheet  of  the  device. 
An  example  of  switching  loss  curve  is  shown  in  Figure  45. 
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Switching  Loss  vs.  Collector  Current 
Edc=300V,Vcc=1 5V,Tj=25°C 


0  50  100  150  200  250  300 

Collector  current :  Ic  (A) 


Figure  45  -  Switching  Loss  vs.  Collector  Current.  The  Total  Switching  Loss  for  one 
Switching  Operation  is  the  Sum  of  the  Three  Shown  in  the  Figure 

We  can  estimate  the  switching  loss  using  the  curve.  For  example,  if  the  load  current  is  100A, 
and  switching  frequency  is  1000Hz,  then  the  switching  loss  for  one  IGBT  is  5  watts.  The  total 
loss  in  the  inverter  is  about  30  watts. 

4.4.2  Power  Loss  on  Unloading  Resistor 

When  the  motor  decelerates,  it  operates  like  a  generator,  and  the  kinetic  energy  is  transformed  to 
electricity.  This  regenerative  energy  will  flow  back  to  the  EU,  and  need  to  be  stored  or 
dissipated. 

To  deal  with  the  regenerative  power,  a  typical  circuit  shown  in  Figure  46  is  often  used.  When 
regenerative  power  is  produced,  current  will  flow  back  from  inverter  and  charge  the  capacitor. 
The  switch  SW1  (a  MOSFET  or  IGBT)  is  off.  If  the  capacitor  can  handle  the  amount  of  energy, 
there  will  be  no  power  dissipation,  and  the  regenerative  energy  will  be  sent  to  the  motor  when 
the  inverter  draws  power  from  DC  bus  to  the  motor. 
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Figure  46  -  Regenerative  Energy  Storage  and  Dissipation  Circuit 

The  capacitor  absorbs  the  regenerative  energy  by  being  charged  by  regenerative  current.  While 
the  capacitor  is  being  charged,  the  DC  bus  voltage  will  increase  due  to  the  capacitor’s  voltage 
equation: 


W  =  r„+£i,'c(<)<A  (74) 

where  Vco  is  the  DC  bus  initial  voltage,  ic  is  the  charging  current,  and  Vc  is  capacitor  voltage. 

Since  the  DC  bus  voltage  has  its  maximum  allowed  value,  the  capacitor’s  ability  to  store  the 
regenerative  energy  is  limited.  If  the  regenerative  energy  exceeds  the  capacitor’s  limit,  SW1  will 
turn  on  and  excessive  power  will  be  dissipated  on  the  unloading  resistor,  producing  power  loss 
and  heat. 

The  regenerative  power  can  be  obtained  by  calculating  the  total  power  input  to  the  motor.  From 
the  EMA  simulation  data  or  experiment  data,  we  have  the  voltage  and  current  values  of  the 
motor.  The  instantaneous  power  consumed  by  the  motor  could  be  calculated: 

PMotor  (0  =  |  Ud  OK  (0  +  h  OK  0)]  (75) 

where  id  and  ud  are  direct  current  and  voltage,  i  and  u  are  quadrature  current  and  voltage. 

When  the  power  input  is  negative,  that  amount  of  regenerative  power  is  generated.  It  is  highly 
dependent  on  the  mission  profile.  However  the  amount  of  regenerative  power  that  is  dissipated 
and  lost  in  the  resistor  is  highly  dependent  on  how  the  circuit  is  designed  and  how  much  energy 
can  be  stored  in  the  capacitor  (e.g.,  the  value  of  the  capacitance  and  the  DC  bus  voltage’s  upper 

ym  )• 

4.4.3  Power  Loss  on  Control  Circuits 

The  control  circuits  process  the  feedback  signals  and  generate  gate  control  signals  to  control  the 
inverter.  The  main  part  is  a  microcontroller  and  its  auxiliary  circuits  and  power  supply.  All  those 
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components  work  in  a  low  voltage  level.  The  power  supply  to  the  microcontroller  is  3.3V  or  5V, 
and  the  typical  power  consumption  is  about  1  to  2  W  depending  on  the  load  of  computation  [24], 

The  working  current  of  those  components  are  at  milliamperes,  and  the  total  power  consumption 
of  the  circuit  is  almost  negligible  compared  to  the  power  loss  on  IGBTs  and  the  unloading 
resistor.  Therefore  when  modeling  the  power  loss  on  the  control  and  drive  board,  the  control 
circuit  part  will  not  be  considered. 

4.5  Modeling  of  the  Power  Loss  on  the  EU 

After  the  analysis  of  the  contributing  power  losses  on  the  EU,  a  model  can  be  built  to  calculate 
the  total  power  loss.  This  total  loss  will  include  the  losses  from  both  contributors:  the  IGBTs  and 
the  unloading  resistor. 

The  power  loss  on  IGBTs  is  calculated  using  the  stator  current,  which  can  be  obtained  from 
motor  simulation  model  [25],  or  from  measurements.  Eq.  70  and  Eq.  72  are  used  to  calculate  the 
power  loss.  Parameters  such  as  RON  ,  Eon ,  Eoff  and  Err  can  be  obtained  from  datasheet  of  the 

device.  If  such  data  is  not  available,  we  can  use  the  values  of  a  device  with  the  same  power 
rating  as  an  approximation. 

The  regenerative  power  is  calculated  using  terminal  voltages  and  currents  of  the  stator,  either 
from  simulation  data  or  measurements.  Follow  Eq.  74  we  can  calculate  the  instantaneous  power 
consumed  by  the  motor,  and  when  it’s  negative  power,  there  is  regenerative  power  flows  through 
the  motor  back  to  the  DC  bus. 

The  DC  bus  capacitor  voltage  is  monitored.  When  the  capacitor  is  being  charged  by  regenerative 
current,  its  voltage  follows  Eq.  73.  When  the  inverter  draws  power  from  the  DC  bus,  the 
capacitor  is  discharged  first  before  drawing  power  from  the  rectifier  and  the  power  system.  It 
follows  the  same  equations  above,  only  the  sign  of  the  current  will  be  negative.  When  the  voltage 
on  the  capacitor  drops  to  the  initial  DC  bus  voltage,  the  EMA  will  be  powered  by  the  rectifier, 
and  the  capacitor  voltage  will  be  locked  at  the  rectifier’s  output  voltage.  To  sum  up,  the  capacitor 
voltage  can  be  written  as: 

if  v„<vc  <  vm 

if  ic  <  0  and  Vc  <  Vco  (76) 

if  ic  >  0  and  Vc  >  Vmax 

When  the  voltage  on  the  capacitor  reaches  the  voltage  limit  (the  third  situation  in  the  above 
equation),  regenerative  current  will  go  through  the  unloading  resistor  and  there  will  be  power 
loss  from  regenerative  power.  The  power  loss  on  the  unloading  resistor  will  be 

Pur  =i2  Run  (77) 

where  i  is  the  regenerative  current  and  Rm  is  the  unloading  resistance. 
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4.6  Thermal  modeling  of  the  EU 


A  Finite  Element  Analysis  (FEA)  based  lumped  node  network  and  its  simulation  of  a  mission 
profile.  This  model  is  based  on  a  detailed  FEA  model  to  locate  the  hot  spots,  to  determine  the 
network  parameters  and  to  verify  its  effectiveness.  The  model  can  also  deal  with  the  nonlinear 
behavior  of  the  EU  introduced  temperature-dependent  magnetic  properties.  This  model  can  be 
incorporated  into  lumped  node  magnetic  and  electric  model  to  develop  a  full  multi-physics, 
multi-scale  simulation  engine.  This  engine  can  accurately  analyze  the  complete  EU  in  a 
systematic  scale  and  whole-mission  duration. 

To  develop  a  lumped  node  thermal  network  to  represent  the  temperature  of  every  solid  part  of 
the  EU.  The  thermal  resistances  and  capacitances  between  the  nodes  can  be  treated  as  electrical 
resistances  and  capacitances.  Hence  the  temperature  of  every  node  can  be  solved  as  the  voltage 
in  this  equivalent  network. 

A  thermal  model  was  present  to  show  the  heat  load  of  the  EU.  The  model  detail  is  based  on 
Danaher  control  unit  (SERVOSTAR  S610).  Figure  41  shows  the  control  unit.  There  are  six 
IGBTs  under  the  heat  spreader  marked  by  “EP”.  The  green  rectangular  show  on  Figure  41  is  the 
regenerative  resistor.  Figure  47  is  the  data  sheet  provided  by  manufacture. 


Model 

S610 

Output  Continuous 

Current  Per  Phase 

{RMS/phase) 

10 

Output  Peak  Current 

Per  Phase  (5  sec) 

20 

Rated  Input  Power 
(KVA)  @  480  VAC 

8.1 

Internal  Power 

90 

Dissipation  Watts 

AC  Input  Line  Voltage 

208-480 

{3  phase) 

VAC 

Continuous  Regen  Power 

Internal  Watts 

200 

External  Watts 

1500 

Figure  47  -  SERVOSTAR  S610  Data  Sheet 

The  model  is  generated  in  FLOTHERM  in  order  to  simulate  the  actual  situation  of  different 
materials,  changing  cross-sections,  and  heat  loads.  In  the  thermal  model,  heat  is  mostly 
generated  by  the  IGBTs,  and  the  model  only  needs  to  simulate  heat  flow  from  IGBTs,  not  from 
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the  rest  of  the  board.  Heat  generated  at  the  IGBTs  is  conducted  to  the  heat  sink  through  the  heat 
spreader  (EP).  A  fan  provides  forced  air  convection  to  remove  heat  from  the  heat  sink. 
According  to  manufacturer  data  sheet,  the  control  board  internal  power  dissipation  is  90  W  at 
480  VAC  power  input.  The  following  assumptions  are  made.  (1)  The  board  is  running  at  half 
power  input  with  240  VAC  and  45  W  power  dissipation;  (2)  80  percent  of  power  goes  to  six 
IGBTs;  (3)  a  Sanyo  Denki  60  x  15  mm,  12  V  axial  fan  provides  forced  air  through  the  heat  sink. 
Figure  48  shows  the  configuration  of  simulation  of  control  unit,  and  the  air  flow  speed  inside  of 
control  unit  (the  top  cover  was  removed  to  show  inside  of  the  control  board). 


FL0THERM 


Figure  48  -  Control  Box  Configurations 


FL0THERM 


FL0THERM 


Figure  49  -  Simulation  Results 

Shown  in  Figure  49,  the  IGBT  temperature  is  about  40°C.  A  commercially  available  fan  is 
adequate  to  manage  the  heat  generated  by  the  control  board. 
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In  order  to  monitor  the  temperature  on  heat  sink,  the  “modified”  simulation  has  been  carried  out. 
There  are  four  temperature  monitor  points  located  on  the  heat  sink  as  shown  on  Figure  50. 


Figure  50  -  Thermocouple  Locations 

AFRL  provided  the  heat  load  data  that  form  the  IGBTs  and  the  regenerative  power.  These  are 
very  high  frequency  data  for  600  seconds  (0.00282s/data).  To  avoid  having  to  use  an  extremely 
small  time  step  for  heat  transfer  calculations,  the  data  is  averaged  by  integrating  over  every  5.6 
seconds.  Since  the  time  scale  for  heat  transfer  is  large,  this  approach  is  acceptable  and  saves 
computation  time. 


Figure  51  -  Regenerative  Power  Load  History 
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Time  (s) 


Figure  52  -  IGBT  Loss  Load  History 

Figure  51  shows  the  time  step  of  regenerative  power  load  history.  Figure  52  shows  the  time  step 
of  IGBT  loss  load  history.  Figure  53  shows  the  Danaher  control  board  simulation  results. 


Figure  53  -  Power  Loss  &  Regenerative  Resistor  Temperature  vs.  Time 


The  simulation  was  carried  out  for  the  entire  600  seconds  of  data  provided  by  AFRL.  The 
temperature  rise  is  very  modest  because  of  the  low  level  of  heat  from  the  IGBTs  and  the  low 
regenerative  power  (less  than  10  W  at  all  time). 
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0.000  0.100(m) 

=  i  Regen-resistor 

Figure  54  -  3-D  Model  of  Control  Unit 

Figure  54  is  a  3-D  model  of  the  Danaher  control  unit.  The  distribution  of  the  lumped  node  is 
listed  in  Table  11,  and  the  lumped  node  network  is  showed  in  Figure  55. 


Table  11  -  Lumped  Node  List 


Node  name 

Material 

Node  number 

IGBT  core 

Silicon 

1 

IGBT  base 

Aluminum  alloy 

2 

Heat  spreader 

Aluminum  alloy 

3 

Heat  sink 

Aluminum  alloy 

4 
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IGBT  base 
2  R2 

— JVVV 

C2  HZ 


1.  IGBT  Core 

2.  IGBT  Base 

3.  Heat  Spreader 

4.  HeatSink 

5.  Regen  resistor 


IS(t) 

0- 


IGBT  core 
1  R1 
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Heat  spreader  Heatsink  Environment  temperature 


Figure  55  -  Lumped  Node  Network 

A  total  thermal  load  give  in  Figure  56  with  a  heat  distribution  36  W/m2  is  assumed.  Figure  57 
show  the  transient  temperature  response  inside  the  EU. 


0  200  400  600  300  1000 

Time  (s) 


Figure  56  -  Thermal  Load 
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Temperature  (degC) 
>  37.94 
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24.485 
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Figure  57  -  One  Pluse  Result 

The  EU  is  assumed  to  have  only  one  thermal  path  with  the  other  parts  and  environment.  The 
thermal  resistnace,  R4,  represents  the  convenction  thermal  resistance  between  heat  sink  and  the 
air  around  it.  Table  12  shows  that  the  thermal  resistator  value  R  and  capciator  value  C  can  be 
calculated  from  steady  state  simulation.  Table  13  shows  that  the  value  of  thermal  resistor  and 
capaciator. 


Table  12  -  Thermal  Resistator  R  and  Capciator  C  Unit 


Thermal  Inputs 

Explanation 

Units 

Source 

Rth 

thermal  resistance 

°C/W 

FEM 

cth 

thermal  capacitance 

J/°C 

FEM 

Tevn 

Environment  temperature 

°C 

Given 

Is 

IGBT  core  loss 

w 

Electrical  Simulation 

Table  13  -  Thermal  Resistator  R  and  Capciator  C  Value 


Model  parameters 

R 

Unit  (°C/W) 

C 

Unit  (J/°C) 

R1 

0.032 

Cl 

1.36 

R2 

0.039 

C2 

10.36 

R3 

0.043 

C3 

48.13 

R4 

0.392 

C4 

860 

R5 

0.392 

C5 

860 
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After  all  the  resistor  and  capacitor  values  of  the  lumped-node  network  model  of  the  EU  are 
know,  the  value  of  the  resistors  and  capaciators  can  be  used  to  simluate  the  temperature  response 
of  the  EU.  A  set  of  equations  can  be  written  as  following: 


(78) 

T\—T2  ^  dT2  T2  —T3  _ 

Rl  2  dt  R2 

(79) 

T2~T1  C  dT3  T3~T4  =  q 

R2  3  dt  R} 

(80) 

T2~Ta  £  dT4  Tevn  —  T4  _ o 
i?3  4  dt  R4 

(81) 

\Tr+<i-QM 

(82) 

Those  equations  set  can  be  solved  by  standard  numerical  methods.  In  Figure  58  and  Figure  59 
are  showing  the  simulation  result  in  the  lumped  node  model  that  compares  with  FEA  result  under 
the  same  initial  condition  and  boundary  conditions. 


Figure  58  -  Result  with  36W  Heat  Load 
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Figure  59  -  Result  with  72W  Heat  Load 

The  result  show  excellent  transient  response  for  the  step  heat  load  given  in  Figure  57.  The 
agreement  between  FEA  results  and  lumped-node  model  proves  that  this  model  is  an  effective 
alternative  to  computational  intensive  FEA  results. 

4.7  An  EU  Example  and  Results 

The  simulation  model  is  built  in  MATLAB.  It  takes  the  electrical  data  and  outputs  the  total 
power  loss  on  the  EU.  The  simulation  takes  a  325  seconds  long  mission  profile.  The  stroke  and 
load  force  profiles  of  the  mission  are  is  shown  in  Figures  60  and  61.  This  mission  profile  is 
highly  dynamic. 
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Figure  60  -  Stroke  of  the  Mission  Profile 


Force 


Figure  61  -  Load  Force  of  the  Mission  Profile 

The  power  loss  in  the  EU  is  calculated  and  plotted  in  Figures  62  and  63.  Because  the  profde  is 
highly  transient,  the  plot  consists  of  lots  of  high  spikes,  which  makes  it  hard  to  see  the  details  of 
the  power  loss.  Further  processing  is  done  to  the  data  to  get  an  average  over  every  second,  so 
that  the  magnitude  of  power  loss  can  be  better  viewed.  Figure  62  is  the  power  loss  on  the 
IGBTs.  Figure  63  is  a  smoothed  version  of  Figure  62.  For  the  regenerative  power  loss,  we  pick 
the  DC  bus  capacitor  to  be  0.1 4F,  and  DC  bus  voltage  limited  to  have  a  340  V  maximum.  Under 
this  setting,  the  capacitor  is  fully  capable  of  buffering  the  regenerative  power,  and  no  power  is 
dissipated  on  the  unloading  resistor  in  this  simulation. 
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Power  loss  on  IGBTs 


Figure  62  -  Power  Loss  on  IGBTs 


Figure  63  -  Power  Loss  on  IGBTs  after  Smoothing 

A  thermal  model  was  developed  to  show  the  temperature  profile  of  the  control  unit.  There  are  six 
IGBTs  under  the  heat  spreader  marked  by  “EP”  on  Figure  41.  Also,  the  large  rectangular  shape 
object  next  to  the  heat  spreader  is  the  unloading  resistor. 

The  model  is  generated  in  FLOTHERM  in  order  to  simulate  the  actual  situation  of  different 
materials,  changing  driver  unit  configurations,  and  heat  loads.  In  the  thermal  model,  heat  is 
mostly  generated  by  the  IGBTs,  and  the  model  only  needs  to  simulate  heat  from  the  IGBTs  and 
the  unloading  resistor,  not  from  the  rest  of  the  board.  Heat  generated  at  the  IGBTs  is  conducted 
to  the  heat  sink  through  the  heat  spreader  (EP).  A  fan  provides  forced  air  convection  to  remove 
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heat  from  the  heat  sink.  In  the  model,  the  assumption  is  made  that  a  Sanyo  Denki  60  x  15  mm, 

12  V  axial  fan  provides  forced  air  through  the  heat  sink. 

The  heat  load  data  for  the  IGBTs  is  provided  in  Figure  62.  These  is  very  high  frequency  data  for 
325  seconds  (0.00282s/datum).  To  avoid  having  to  use  an  extremely  small  time  step  for  heat 
transfer  calculations,  the  one-second  averaged  data  shown  in  Figure  63  is  used. 

A  thermal  simulation  is  done  using  the  power  loss  data  from  our  model,  and  the  temperature  of 
the  IGBTs  is  calculated.  An  initial  temperature  of  28°C  and  an  ambient  temperature  of  26°C 
have  been  applied  to  the  model  to  simulate  the  control  unit.  A  comparison  between  the 
simulation  and  measurement  of  the  IGBTs  is  shown  in  Figure  64.  It  shows  that  the  simulation 
agrees  well  with  the  measured  temperature  on  the  device.  The  IGBTs’  temperature  remains 
approximately  2  °C  higher  than  the  ambient  temperature  which  is  sufficient  to  remove  the  heat 
from  the  IGBT  package.  A  comparison  for  the  unloading  resistor  is  shown  in  Figure  65.  In  this 
mission  profile,  the  unloading  resistor  is  barely  used,  the  capacitor  is  sufficient  to  handle  the 
regenerative  power,  and  very  little  power  is  dissipated  on  the  resistor. 


IGBT  Temperature 


Time  (s) 

Figure  64  -  Temperature  on  IGBTs.  Hardware  Test  Data  is  compared  with  Simulation 

Data  Based  on  the  Power  Loss  Model 
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Figure  65  -  Temperature  on  Resistors.  We  Pick  the  DC  Bus  Capacitor  to  be  0.14F,  and 
DC  Bus  Voltage  Limited  to  340V  Maximum  for  the  Regenerative  Power  Loss  calculation.  It 
shows  that  the  Unloading  Resistor  is  barely  used  and  the  Capacitor  is  Fully  Capable  of 
Absorbing  Regenerative  Power  for  this  Mission  Profile 

The  power  loss  on  the  EU  board  is  analyzed  and  modeled  using  MATLAB.  Most  of  the  power 
loss  on  the  EU  comes  from  the  IGBTs  in  the  power  inverter,  and  the  regenerative  power 
dissipation  resistor.  The  switching  loss  and  on-state  loss  are  considered  for  the  IGBTs. 

Datasheets  are  used  to  obtain  the  input  parameters  of  the  model.  A  simulation  is  done  in 
MATLAB,  and  the  power  loss  data  is  plotted.  Thermal  analysis  is  done  and  the  results  are 
compared  with  the  measured  temperatures.  The  agreement  validates  the  power  loss  model  for  the 
IGBTs  and  the  unloading  resistor. 
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5.  CONCLUSION 


There  is  a  strong  motivation  to  move  towards  all-electric  aircraft.  However,  the  significant 
issues  to  full  realization  of  all-electric  aircraft  are  the  thermal  management  of  heat  loads  and 
power  management.  We  have  successfully  developed  and  implemented  an  EMA  model  with  a 
control  algorithm,  a  PWM  component,  an  electromechanical  component  with  nonlinear 
parameters,  and  a  thermal  component.  The  model  takes  a  desired  position  profile,  a  load  force 
profile,  and  an  environmental  temperature  profile  as  inputs  and  predicts  the  thermal  behavior  of 
the  motor  and  the  transient  power  demands  of  the  EMA.  It  was  shown  that  the  heat  generation 
profiles  are  indeed  highly  dynamic,  and  control  plays  a  powerful  role  in  how  much  excess  heat  is 
generated  and  how  much  the  power  demand  spikes. 

The  lumped-node  thermal  network  is  a  feasible  approach  to  simulate  the  EMA  thermal 
characteristics.  A  thermal  network  is  based  on  tuning  the  thermal  resistance  and  thermal 
capacitance  values  with  the  FEM  results.  The  lumped-node  model  has  been  successfully  applied 
to  step  and  transient  heat  loads  with  cooling  by  both  convection  and  radiation.  From  the  cases 
studied,  it  could  be  concluded  that  the  major  barrier  to  heat  removal  from  the  motor  is  external  to 
the  motor.  Strong  forced  convection  is  needed.  For  situations  where  phase  change  materials 
may  be  needed  to  store  the  thermal  energy  to  avoid  the  motor  reaching  excessively  high 
temperature,  a  PCM  module  has  been  incorporated  successfully  into  the  lumped  node  model. 

The  power  loss  in  the  electronic  unit  (EU)  of  an  EMA  system  has  been  simulated  also.  The 
control  and  drive  unit  (EU)  processes  the  input  data  and  feedback  from  sensors,  and  deliver  input 
voltage  and  current  to  the  motor  through  an  insulated-gate  bipolar  transistor  (IGBT)  module. 
Most  of  the  power  loss  on  EU  comes  from  the  IGBTs  in  the  power  inverter,  and  the  regenerative 
power  dissipation  resistor.  The  switching  loss  and  on-state  loss  are  considered  for  IGBTs. 

The  model  was  also  used  to  compare  its  temperature  predictions  with  temperature  measurements 
of  the  motor  and  power  drive  provided  by  the  Air  Force  Research  Laboratory.  Though  the 
agreement  was  good,  due  to  the  low  heat  load  associated  with  the  test  cases,  the  model  could  not 
be  validated  conclusively. 
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APPENDIX  A:  HOW  TO  GET  PARAMETERS  FOR  AN  EMA 

SIMULATION 

A 1  What  Parameters  Are  needed 

A  simulation  model  of  an  electromechanical  actuator  (EMA)  requires  accurate  parameters  in 
order  to  give  accurate  simulation  results.  The  primary  parameters  that  are  required  for  the 
simulation  model  to  run  are  listed  here: 

•  Gearing  ratio 

•  Number  of  poles 

•  Power  train  friction 

•  Shaft  mass 

•  Moment  of  inertia 

•  Phase  resistance 

•  PM  flux  linkage 

•  Direct  inductance 

•  Quadrature  inductance 

•  Thermal  capacitance 

•  Thermal  resistance 

These  parameters  are  necessary  for  the  EMA  simulation  model.  In  addition  to  these  parameters, 
it  is  helpful  to  know  the  limits  of  the  motor  and  actuator  assembly.  These  limits  include 

•  No-load  speed 

•  Stall  force 

•  Stroke  range 

These  three  have  less  to  do  with  how  the  EMA  itself  behaves  and  more  with  how  it  should  be 
controlled.  For  example,  there  is  no  direct  restriction  on  how  much  force  an  EMA  must  face 
before  it  stalls,  but  a  force  greater  than  a  certain  value  would  require  more  current  than  the 
windings  of  the  motor  could  safely  handle  without  burning  out.  The  motor  controller  prevents 
this  from  happening,  so  we  get  the  effect  of  a  limit  to  the  force  we  can  see  from  the  EMA. 
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A2  Where  Parameters  Can  Be  Found 

Some  of  these  parameters  are  easy  to  determine  and  others  require  careful  analysis.  There  are 
two  approaches  to  determining  these  parameters.  One  is  to  take  a  physical  EMA  and  run  various 
experiments  on  it  measuring  the  parameters  one-by-one.  This  is  often  a  destructive  process  since 
it  ultimately  requires  completely  taking  apart  the  machine  in  order  to  see  its  internal  construction. 
The  other  approach  is  to  design  and  model  a  motor  that  is  as  similar  to  the  target  motor  as 
possible  and  use  the  parameters  from  the  model  for  the  full  EMA  simulation.  This  is  non¬ 
destructive,  but  has  the  potential  for  somewhat  higher  uncertainties.  There  are  some  parameters, 
however,  that  can  be  measured  without  taking  any  components  of  the  EMA  apart. 

A2.1  What  Can  Be  Found  by  Measurement 

The  following  sections  explain  how  the  various  parameters  of  an  EMA  can  be  found 
experimentally. 

A2.1.1  No-load  speed 

The  no-load  speed  can  be  found  by  sending  position  commands  to  the  actuator  to  move  back  and 
forth.  This  is  done  without  any  load  attached  to  the  EMA.  Taking  the  minimum  of  the  speeds  in 
the  two  directions  gives  the  no-load  speed. 

A2.1.2  Stall  force 

The  stall  force  is  measured  while  the  actuator  is  connected  to  a  hydraulic  press  or  other  active 
load.  The  load  force  from  the  hydraulic  press  is  increased  while  the  actuator  is  commanded  to 
hold  its  position  until  the  actuator  begins  to  move.  When  the  actuator  begins  to  move,  it  is  no 
longer  able  to  sustain  the  force  of  the  load,  and  this  is  the  stall  force. 

A2.1.3  Stroke  range 

The  stroke  range  is  perhaps  one  of  the  easiest  parameters  to  determine.  It  can  come  from  either 
design  specifications  or  from  simply  extending  the  actuator  as  far  as  it  can  reasonably  go. 

A2. 1 .4  Gearing  ratio 

Depending  on  the  documentation  from  the  manufacturer,  the  gearing  ratio  can  be  found  either 
from  the  documentation  of  the  actuator  or  by  taking  the  gear  train  apart,  counting  the  teeth  on 
each  gear  and  measuring  the  linear  displacement  of  the  actuator  rod  for  one  full  revolution  of  the 
ball  screw.  Although  this  does  involve  some  disassembly  of  the  EMA,  it  is  generally  a  non¬ 
destructive  process. 

A2. 1 .5  Number  of  poles 

The  documentation  of  the  motor  in  the  EMA  should  state  the  number  of  poles.  If  not,  it  can  be 
found  either  by  taking  the  motor  apart  or  by  scanning  it  with  a  high-intensity  X-ray  machine. 
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A2. 1 .6  Power  train  friction 


One  method  of  finding  the  power  train  friction  involves  using  a  linear  load  such  as  a  hydraulic 
press  coupled  to  the  actuator.  The  EMA  is  driven  forward  and  backward  while  the  EMA  is 
unpowered.  Using  the  force  and  stroke  data,  the  static  friction  can  be  obtained.  First,  the  mean 
forces  required  to  push  or  pull  the  EMA  rod  are  found.  Depending  on  the  orientation  of  the  EMA, 
the  magnitudes  of  these  forces  will  not  be  equal  due  to  a  bias.  The  bias  exists  because  of  the 
weight  of  the  rod  of  the  EMA.  The  difference  between  this  bias  and  the  push  and  pull  forces  is 
the  friction  force.  As  a  dry  (or  contact)  friction,  it  is  mostly  constant  over  a  wide  range  of  speeds. 

A2.1.7  Shaft  mass 

The  mass  of  the  whole  shaft  of  the  EMA  can  be  derived  from  the  friction  analysis.  The  bias  in 
the  force  required  to  move  the  EMA  rod  is  due  to  the  weight  of  the  rod. 

A2.1.8  Moment  of  inertia 

tia  of  the  motor’s  rotor,  the  gears,  and  the  ball  screw  can  be  found  by  accelerating  the 
actuator  under  impulse  from  the  hydraulic  press.  The  curvature  of  the  velocity  curve  over  time  is 
a  function  of  the  total  inertia  of  the  drive  train.  This  inertia  less  the  mass  of  the  shaft  is  the 
moment  of  inertia. 

Determining  these  two  inertias  (shaft  mass  and  moment  of  inertia)  separately  is  important  in 
order  to  get  an  accurate  force  of  gravity  and  net  force  in  the  simulations. 

A2.1.9  Phase  resistance 

Of  all  the  electrical  parameters,  the  phase  resistance  is  the  easiest  to  obtain.  It  is  measured  by 
connecting  leads  to  two  of  the  three  lines  into  the  motor  and  measuring  the  open-circuit  line-to- 
line  resistance.  This  value  is  then  divided  by  two  if  the  motor  is  a  Wye-connected  motor.  The 
resulting  value  is  the  phase  resistance  of  the  EMA.  Any  pair  of  two  lines  should  measure  the 
same  resistance.  If  they  do  not,  there  is  likely  a  short  in  the  motor. 

The  remaining  parameters  (PM  flux  linkage,  direct  and  quadrature  inductances,  and  thermal 
resistances  and  capacitances)  are  very  difficult  to  measure  without  destructively  taking  the 
electric  motor  apart. 

A2.2  What  Can  Be  Found  by  Modeling 

The  following  parameters  are  much  more  readily  obtained  by  modeling  than  by  experimental 
measurement: 

•  PM  flux  linkage 

•  Direct  inductance 

•  Quadrature  inductance 
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•  Thermal  capacitance 

•  Thermal  resistance 

Because  of  the  difficulty  of  measuring  these  parameters  directly,  developing  a  model  and 
calculating  these  parameters  from  the  model  can  be  very  helpful.  The  procedure  for  developing 
such  a  model  is  detailed  in  the  following  section. 

A3  How  to  Model  a  Motor  for  Parameter  Estimation 

There  are  very  standard  designs  for  electric  motors.  A  given  speed  and  power  requirement  can  go 
a  long  way  in  deciding  the  specific  design  of  a  motor.  This  fact  makes  it  reasonable  to  using 
good  design  process  to  get  a  reasonable  estimate  of  the  five  remaining  parameters  mentioned  in 
the  previous  section: 

•  PM  flux  linkage 

•  Direct  inductance 

•  Quadrature  inductance 

•  Thermal  capacitance 

•  Thermal  resistance 

There  are  four  basic  steps  to  getting  these  parameters  from  a  model.  The  steps  are 

•  Design  an  electric  motor  to  specification 

•  Refine  the  design 

•  Model  the  motor  with  finite  element  method  (FEM) 

•  Extract  the  reduced-order  parameters 
The  following  sections  will  detail  these  steps. 

A3.1  Design  an  electric  motor  to  specification 

To  design  an  electric  motor  some  basic  specifications  are  required,  such  as  the  input  voltage  the 
motor  will  see  or  the  number  of  poles  the  motor  should  have.  Included  with  the  software  package 
to  AFRL  is  a  script  called  “pmdesigner.m”.  This  program  takes  the  user  step-by-step  through  the 
process  of  designing  a  permanent  magnet  (PM)  electric  motor.  A  sample  of  the  questions  asked 
by  pmdesigner  is  shown  here: 

Do  you  want  to  start  a  new  design?  [y,N] :  y 

Enter  the  name  of  the  new  design  txt  file:  myDesign 

Pout  (Output  power  expected)  bounded  by  0  and  Inf  [W] :  300 
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VT  (Terminal  voltage)  bounded  by  0  and  Inf  [V] :  120 
nm  (Mechanical  speed)  bounded  by  0  and  Inf  [krpm] :  1.6 
DY  (Connection  type:  D  or  Y)  limited  to  one  of  D,  Y:  y 

The  value  given  for  DY  does  not  fit  the  specified  limits 
DY  (Connection  type:  D  or  Y)  limited  to  one  of  D,  Y:  Y 
p  (Number  of  poles)  bounded  by  0  and  Inf:  12 
sp  (Fractional  stator  pitch)  bounded  by  0  and  1:  1 
Ns  (Target  number  of  slots)  bounded  by  0  and  Inf:  33 
zeta_m  (Mechanical  skew)  bounded  by  0  and  Inf:  0 
Br  (Remanence)  bounded  by  0  and  Inf  [G] :  1.4 
He  (Coercivity)  bounded  by  0  and  Inf  [Oe] :  9.31e5 
eff  (Efficiency  (ratio,  not  percentage) )  bounded  by  0  and  1 
0.95 

pf  (Power  factor  (ratio))  bounded  by  0  and  1:  1 
pfSgn  (Power  factor  sign  (1  or  -1))  limited  to  one  of  -1,  1 

1 

osf  (Overspeed  factor  (ratio))  bounded  by  0  and  1:  0.2 
emb  (Magnet  embrace)  bounded  by  0  and  1:  0.95 
slotStyle  (Slot  style  (box:  1,  tight:  2,  wedge:  3))  limited 
to  one  of  1,  2,  3:  2 

useLining  (Use  lining  (1)  or  not  (0))  limited  to  one  of  0, 

1:  1 

coolingMethod  (Cooling  method  (air:  1,  water:  2))  limited 
to  one  of  1 ,  2 :  1 

The  details  of  this  particular  setup  will  be  saved  to  a  file  called  “myDesign,”  as  specified  above. 
A  design  of  an  electric  motor  using  several  carefully  chosen  design  equations  is  then  output  by 
pmdesigner.  The  results  for  this  particular  setup  are  shown  here: 


Pout  = 

300 

w 

Output  power 

Pfw  = 

10 

w 

Friction  loss 

Pair  = 

1.0089 

w 

Air  friction  loss 

Pcore  = 

7 . 8  92  6e-0  6 

w 

Core  loss 

VT  = 

120 

V 

Rated  voltage 

Vph  = 

69.282 

V 

Phase  voltage 

P  = 

12 

Poles 

Np  = 

2.75 

Slots  per  pole 

Nm  = 

3 

Slots  spanned  by  a  coil 

Ns  = 

33 

Slots  on  stator 

nm  = 

1600 

rpm 

Speed 

fm  = 

26.667 

Hz 

mechanical  frequency 
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fe 

= 

160 

Hz 

electrical  frequency 

D 

= 

59.02 

mm 

Stator  bore  diameter 

DO 

= 

43426 

mm 

Stator  core  diameter 

g 

= 

28 . 809 

mm 

Actual  air  gap 

ge 

= 

29.51 

mm 

Effective  air  gap 

Dr 

= 

1 .4016 

mm 

Actual  rotor  diameter 

vr 

= 

23 . 114 

f t/min 

Actual  rotor  tip  speed 

0 .066039 

o 

o 

of  maximum  allowable 

1 

= 

59.02 

mm 

Rotor  length 

Nc 

= 

127741 

Turns  per  coil 

Ne 

= 

2 . 5513e+0  6 

Effective  turns  per  phase 

ef f  N 

= 

181 . 57 

o 

o 

Turns  efficiency 

lambda  s 

= 

5.6187 

mm 

Slot  arc  spacing 

bs 

= 

4.45 

mm 

Slot  width 

ts 

= 

1 .1687 

mm 

Tooth  width 

bsO 

= 

1.4833 

mm 

Mouth  Width 

dsl 

= 

21677 

mm 

Slot  depth  (not  including 

wedge) 

dsOb 

= 

1.016 

mm 

Trapezoidal  section  of 

wedge 

depth 

dsOa 

= 

1 . 016 

mm 

Mouth  section  of  wedge 

depth 

sw 

= 

2.032 

mm 

Slot  wedge  depth 

si 

= 

0 . 11477 

mm 

Slot  lining  thickness 

Ss 

= 

96468 

mmA2 

Total  slot  area 

Sc 

= 

91486 

mmA2 

Total  slot  area  minus 

lining 

and  wedge 

a 

= 

1.0545 

Slot  inflation 

Sa 

= 

0 . 18993 

mmA2 

Conductor  cross-sectional 

area 

Da 

= 

0.49176 

mm 

Conductor  diameter 

awg 

= 

24 . 324 

Wire  gauge 

Vol 

= 

8.7414e+10  mmA3 

Total  volume  of  motor 

Malformed  stator  slots:  dsl  is  too  large.  Try  reducing  the 
number  of  slots  on  the  stator. 

Notice  that  pmdesigner  followed  the  request  for  33  slots  but  indicates  that  this  would  result  in  a 
poor  design.  This  means  that  another  pass  should  be  done  with  a  value  smaller  than  33. 

Once  the  design  is  settled  upon,  the  rest  of  the  output  results  will  be  helpful  for  modeling  the 
motor.  For  the  above  setup,  the  rest  of  the  results  are  shown  here: 
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Converted  to  RMXprt  variables: 
Machine : 

Number  of  poles  12 

Frictional  loss  10  W 


Windage  loss 

1 . 0089  W 

Reference  speed 

1600  rpm 

Stator : 

Outer  diameter 

43426  mm 

Inner  diameter 

59.02  mm 

Length 

59.02  mm 

Number  of  Slots 

33 

Slot  type 

3 

Slot: 


HsO 

1 . 016 

mm 

Hsl 

1 . 016 

mm 

Hs2 

21677 

mm 

BsO 

1 .4833 

mm 

Bsl 

4.45 

mm 

Bs2 

4.45 

mm 

Rs 

0 

mm 

Winding : 

Winding  layers 
Winding  type 
Parallel  branches 
Conductors  per  slot 
Coil  pitch 
Number  of  strands 
Wire  wrap 
Wire  size 
Slot  lining 
Wedge  thickness 
Layer  insulation 
Limited  fill  factor 
Rotor : 

Outer  diameter 
Inner  diameter 
Length 
Pole  type 
Pole : 


whole 


0 


0  . 


1.4016  mm 
0.66203  mm 
59.02  mm 
2 


2 

coiled 

1 

255482 

3 

1 

0.0635 

24.5 

.11477 

2.032 

057384 

0.999 


mm 

awg 

mm 

mm 

mm 
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mm 


Embrace 
Offset 
Magnet  type 
Magnet  thickness 
Shaft :  non-magnetic 


0.95 
0 

SmCo28 
0.038753  mm 


These  specifications  will  then  be  entered  into  a  program  by  Ansys  called  RMxprt.  It  is  a  motor 
designing  program. 

A3. 2  Refine  the  design 

RMxprt  is  a  motor  designing  program  that  will  take  the  rough  design  specifications  developed 
above  and  fill  in  the  details.  It  provides  the  bridge  between  the  simple  table  of  design  outputs  in 
the  previous  section  to  a  full-blown  finite  element  model.  The  steps  to  model  a  motor  using 
RMxprt  are  detailed  below  in  a  succinct  manner  step-by-step. 

A3.2.1  Create  new  motor  project 

A3. 2. 1 . 1  Create  new  project 


1  Ansoft  Maxwell  -  Projectl  -  [Slotsl2polesl0  - 


Pro: 


File 


Edit  View  Project  Tools  Windo 


C)  New 


Ctrl+N 


E?  Open.., 
Close 


Ctrl+O 


Q  Save 
Save  As... 

Save  As  Technology  File 


Ctrl+S 


Page  Setup... 

Figure  A-l  -  Selecting  “New”  from  the  “File”  Menu 

A3.2.1.2  Insert  RMxprt  Design 


Project  Tools  Window  Help 


^  Insert  Maxwell  3D  Design 
H|  Insert  Maxwell  2D  Design 


Insert  RMxprt  Design 


Insert  Documentation  File ... 


Analyze  All 


Project  Variables... 

Datasets... 

Remove  Unused  Definitions . 


Figure  A-2  -  Creating  a  New  RMxprt  Design  from  the  “Project”  Menu 

77 

Approved  for  public  release;  distribution  unlimited. 


A3.2.1.3  Select  Brushless  Permanent-Magnet  DC  Motor 


Select  Machine  ype 


Three-Phase  Induction  Motor 
fijr Single-Phase  Induction  Motor 
iji Three-Phase  Synchronous  Machine 


Brushless  Permanent-Magnet  DC  Motor 


1^1  Adjust-Speed  Synchronous  Machine 
^Permanent-Magnet  DC  Motor 
^Switched  Reluctance  Motor 
.(Si Line-Start  Permanent-Magnet  Synchronous  Motor 
£  Universal  Motor 
ODC  Machine 
Claw-Pole  Alternator 

^Three-Phase  Non-Salient  Synchronous  Machine 

Figure  A-3  -  List  of  Machine  Types 

A3.2.1.4  Save  the  project  and  name  it 
A3. 2.2  Define  machine  properties 

Under  Project  Manager,  select  the  components  listed  below  one-by-one,  each  time  specifying  the 
properties  shown  in  the  Properties  window. 


Project  Manager 


l=Hjp  Project  T 

RMxprt  Design  1  (Brushless  Permanent-Magm 

|  l 

! . ®  Dncuit 

□-Q  Stator 
|  l|  Slot 

! . (JJ)  Winding 

|j-4>  Rotor 

j--CD  Pole 

Shaft 

Analysis 

j--^|  Optometries 
! . HI  Results 


Figure  A-4  -  Project  Manager  Tree 
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|  Properties  *  x 


Name 

Value 

Unit 

Evaluated  Value 

Outer  Diameter 

77.1 

mm 

77.1mm 

Inner  Diameter 

43.5 

mm 

43,5mm 

Length 

53.7 

mm 

5S.7mm 

Stacking  Factor 

0  95 

Steel  Type 

steel 101£| 

Number  of  Slots 

12: 

Slot  Type 

3  J 

Skew  Width 

0 

0 

Figure  A-5  -  Properties  List 


A3.2.2.1  Machine 

•  Machine  Type:  The  machine  type  you  selected  when  inserting  a  new  RMxprt  design, 
Three  Phase  Synchronous  Machine.  It  is  already  specified  by  having  selected 
Brushless  Permanent-Magnet  DC  Motor. 

•  Number  of  Poles:  The  number  of  poles  the  machine  contains.  This  value  is  the  total 
number  of  poles  in  the  stator  (or  the  number  of  pole  pairs  multiplied  by  two). 

•  Frictional  Loss:  the  frictional  power  loss  pfr  (due  to  friction)  measured  at  the 
reference  speed  nr.  The  frictional  power  loss  pf  at  a  speed  n  is  given  by  pf  =  pfr  * 
(n/nr). 

•  Windage  Loss:  The  wind  loss  pwr  (due  to  the  air  resistance)  measured  at  the 
reference  speed  nr.  The  wind  power  loss  pw  at  the  real  speed  n  is  given  by  pw  =  pwr 
*  (n/nr)A3. 

•  Reference  Speed:  The  reference  speed  at  which  the  frictional  power  loss  and  windage 
power  loss  are  measured.  It  is  the  main  speed  the  machine  is  expected  to  run. 

•  Control  Type:  (typical  is  DC) 

•  Circuit  Type:  (typical  is  Y3) 

A3.2.2.2  Circuit 

•  Lead  Angle  of  Trigger:  The  trigger’s  lead  angle  in  electrical  degrees.  A  positive  value 
represents  a  lead  angle,  and  a  negative  value  represents  a  lag  angle.  An  angle  of  0 
means  that  the  average  induced  voltage  in  the  triggered  phase  during  the  Trigger 
Pulse  Width  is  at  a  maximum. 

•  Trigger  Pulse  Width:  The  period  from  the  on-status  to  the  off-status  of  a  transistor,  in 
electrical  degrees. 
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•  Transistor  Drop:  The  voltage  drop  across  one  transistor  when  the  transistor  is  turned 
on. 


•  Diode  Drop:  The  voltage  drop  of  one  diode  in  the  discharge  loop.  If  you  selected  a 
star-type  circuit  (S3  or  S4)  as  the  Circuit  Type,  this  is  the  total  discharge  voltage. 

A3.2.2.3  Stator 

•  Outer  Diameter:  The  outer  diameter  (or  core  diameter)  of  the  stator  core. 

•  Inner  Diameter:  The  inner  diameter  (or  bore  diameter)  of  the  stator  core. 

•  Length:  The  length  of  the  stator  core. 

•  Stacking  Factor:  The  stacking  factor  for  the  stator  core.  It  is  the  ratio  of  metal  (versus 
lamination)  thickness  to  total  thickness.  1  indicate  no  lamination.  0.95  is  common. 

•  Steel  Type:  the  stator  core  material.  This  must  be  of  type  Steel.  Some  good  selections 
are  steel-1010  or  iron. 

•  Number  of  Slots:  The  total  number  of  slots  the  stator  core  contains. 

•  Slot  Type:  The  type  of  the  slots  in  the  stator  core.  Slot  types  1  though  4  are  filled  with 
round  wire.  Slot  types  5  and  6  are  filled  with  rectangular  wire. 

•  Skew  Width:  The  slot  skew  width,  measured  in  slot  number.  This  means,  as  the 
laminations  are  placed  one  on  top  of  the  other,  they  might  not  be  positioned  exactly 
the  same  each  time;  so,  from  lamination  to  lamination,  the  slots  skew,  or  twist,  about 
the  center  of  the  center  of  the  machine.  The  figure  here  shows  a  slot  skew  of  1 . 
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Figure  A-6  -  Model  of  Stator  with  Skewed  Slots 

A3. 2.2.4  Slot 

•  Auto  Design: 

•  Parallel  Tooth: 

•  Tooth  Width: 

•  HsO:  Slot  opening  height 

•  HsOl :  Slot  closed  bridge  height 

•  Hsl :  Slot  wedge  height 

•  Hs2:  Slot  body  height 

•  BsO:  Slot  opening  width 

•  Bs  1 :  Slot  wedge  maximum  width 

•  Bs2:  Slot  body  bottom  width,  0  for  parallel  teeth 

•  Rs:  Slot  body  bottom  fillet  for  slot  types  3  and  4 
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A3.2.2.5  Winding:  Winding 

•  Winding  Layers:  The  number  of  coil  layers  (1  or  2)  in  the  stator  winding.  This  does 
not  depend  on  whether  the  machine  has  lap  winding  or  concentric  winding. 

•  Winding  Type:  The  type  of  the  stator  winding.  To  specify  your  own,  select  Editor, 
right  click  in  the  table  to  the  right,  and  select  Edit  Layout... .  T  means  top  layer 
(closest  to  the  mouth  of  the  slot)  and  B  means  bottom  layer.  For  concentric  winding, 
just  select  Whole. 

•  Parallel  Branches:  The  number  of  parallel  branches  in  one  phase  of  the  stator  winding. 
If  all  the  windings  of  a  phase  are  connected  in  series,  then  the  number  of  parallel 
branches  is  1 . 

•  Conductors  per  Slot:  The  total  number  of  conductors  in  each  stator  slot.  This  value  is 
the  number  of  turns  per  coil  multiplied  by  the  number  of  layers.  Enter  0  to  have 
RMxprt  auto-design  this  value. 

•  Coil  Pitch:  (This  is  enabled  when  Winding  Layers  is  2.)  It  is  the  number  of  teeth 
separating  two  sides  of  one  coil.  For  example,  if  a  coil  starts  in  slot  1  and  ends  in  slot 
2,  it  has  a  coil  pitch  of  one.  For  concentric  winding,  coil  pitch  should  be  1 . 

•  Number  of  Strands:  The  number  of  wires  per  conductor.  Enter  0  to  have  RMxprt 
auto-design  this  value.  For  simple,  small  motors,  this  is  typically  just  1. 

•  Wire  Wrap:  The  thickness  of  the  double-sided  wire  wrap  (electrical  insulation).  Enter 
0  to  automatically  obtain  this  value  from  the  wire  specification  library. 

•  Wire  Size:  Wire  diameter  (for  round  wires)  or  width  and  thickness  (for  rectangular 
wires).  (0  for  auto-design).  You  can  assign  wire  size  of  round  wires  or  rectangle  wires. 
When  the  slot  type,  under  Stator,  is  1  to  4,  round  wires  are  used.  When  the  slot  type 
you  selected  is  5  or  6,  rectangle  wires  are  used. 

A3.2.2.6  Winding:  End/Insulation 

•  Input  Half-turn  Length:  Select  or  clear  this  check  box  to  specify  whether  or  not  you 
want  to  enter  the  half-turn  length  manually. 

•  Half-turn  Length:  The  half-turn  length  of  the  armature  winding.  It  is  available  when 
Input  Half-turn  Length  is  selected. 

•  End  Adjustment:  the  distance  of  one  end  of  a  conductor  extending  vertically  beyond 
the  end  of  the  stator.  It  is  available  when  Input  Half-turn  Length  is  cleared.  Enter  0 
for  auto-design. 

•  Base  Inner  Radius:  The  inner  radius  of  the  base  comer:  how  big  of  an  arc  the  coil 
makes  at  the  end.  Enter  0  for  auto-design. 
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•  Tip  Inner  Diameter:  The  inner  diameter  of  the  coil  tip,  seen  as  Top-End  Inner 
Diameter  above. 

•  End  Clearance:  The  end  clearance  between  two  adjacent  stator  coils. 

•  Slot  Liner:  The  thickness  of  the  slot  liner  insulation. 

•  Wedge  Thickness:  The  thickness  of  the  wedge  insulation. 

•  Layer  Insulation:  The  thickness  of  the  insulation  layer  which  is  enabled  when 
preformed  coils  (for  slot  types  5  or  6)  are  used. 

•  Limited  Fill  Factor:  The  limiting  slot  fill  factor  for  the  wire  design.  This  value  relates 
to  how  much  stuff  is  put  into  the  slot  versus  the  total  space  available  in  the  slot. 


Inner  Diameter 


^Wedge 


Liner 


Coil 


Coil  Wrap 


Layer 

Insulation 


Bottom 

Insulation 


Figure  A-7  -  End  Winding  Characteristics  (Left)  and  Slot  Filling  (Right) 


A3.2.2.7  Rotor 

•  Outer  Diameter:  the  outer  diameter  of  the  rotor  core,  includes  the  magnets. 

•  Inner  Diameter:  the  inner  diameter  of  the  rotor  core;  diameter  of  axle. 

•  Length:  The  length  of  the  rotor  core. 

•  Steel  Type:  This  should  be  a  steel  type,  just  like  with  the  stator. 

•  Stacking  Factor:  The  stacking  factor  for  the  stator  core.  It  is  the  ratio  of  metal  (versus 
lamination)  thickness  to  total  thickness.  1  indicates  no  lamination.  0.95  is  common. 

•  Pole  Type:  This  is  the  arrangement  of  the  magnets  on  the  rotor. 

A3.2.2.8  Pole 


•  Embrace:  This  is  a  ratio  of  the  arc  length  of  a  permanent  magnet  to  the  available  arc 
length  for  each  magnet.  So,  1  means  that  the  magnets  are  all  touching  side-to-side. 
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•  Offset:  The  distance  from  the  center  of  the  rotor  that  the  circle  that  forms  the  outer 
edge  of  the  magnet  is  positioned.  This  is  illustrated  in  this  figure: 


Figure  A-8  -  Diagram  of  Magnet  Arc  Offset  from  Center  of  Axel 

•  Magnet  Type:  The  material  that  the  magnets  are  made  of. 

•  Magnet  Thickness:  The  distance  from  the  inner  arc  to  the  outer  arc  of  the  magnet  at 
its  thickest  point. 

A3.2.2.9  Shaft 


Magnetic  Shaft:  A  check  box  for  whether  the  axle  should  be  made  of  soft  magnetic  material  or  of 
non-magnetic  material.  In  other  words,  should  the  axle  have  a  permeability  equal  to  air  or  should 
it  be  much  greater  than  that  of  air. 


A3. 2. 3  Analyze 

A3.2.3.1  Select  Add  Solution  Setup...  on  Analysis 


Shaft 


& 

■1^1  Optir 


i . [M|  Rest 

m  Definitior 
g]  PFj  Mate 


Paste 

^  Add  Solution  Setup... 
List... 

Analyze  All 


Figure  A-9  -  Selecting  “Add  Solution  Setup...”  from  “Analysis” 
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A3. 2. 3. 2  Enter  the  conditions  under  which  you  wish  the  machine  to  operate 


General  j  Defaults  | 


Operation  Type: 

|  Motor 

"3 

Load  Type: 

Rated  Output  Power: 

F 

[w- 

T  | 

Rated  Voltage: 

1° 

F~ 

T  | 

Rated  Speed: 

|o 

Jrpm 

— 3 

Operating  Temperature: 

|75 

fcei 

Figure  A-10  -  Solution  Setup  Ratings 


•  Rated  Output  Power:  The  output  mechanical  power  developed  at  the  shaft  for  the 
motor. 

•  Rated  Voltage:  The  RMS  line-to-line  voltage. 

•  Rated  Speed:  The  desired  output  speed  of  the  motor  at  the  load  point 

•  Operating  Temperature:  The  temperature  at  which  the  system  functions,  and  select 
the  units.  The  Operating  Temperature  will  affect  all  winding  resistances  and  therefore 
affect  all  ohmic  losses. 


A3.2.3.3  Analyze  the  solution  setup 


I  Copy 
Paste 


Analysis 

Optime 
_  Resuit! 

m  Definitions  Rename 
i-D  Materia  x  pe!ete 
S  Nd 
■  ste 


Properties.. 


Analyze 


Create  £>uick  Report... 


Performance... 
Design  Sheet... 
Jp  Curve... 


Figure  A-ll  -  Selecting  66 Analyze”  from  “Setup  1” 
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A3.2.3.4  Look  at  results.  On  results,  select  Solution  Data... 


Setup  1 
]  Optimetrics 


M  Definitic  Paste 

\k-M  Mai 


Create  RMxprt  Report 
Delete  All  Reports 
Report  Templates 

Output  Variables... 
Update  All  Reports 
Open  All  Reports 

C  reate  Quick  Rep  o  rt. . . 


Solution  Data... 


Figure  A-12  -  Selecting  “Solution  Data...”  from  “Results’ 

This  will  bring  up  a  window  like  this: 

|  Setupl 


Simulation: 


Design  Variation: 


^  |  |  Performance 


7E1 


j  W 


Performance  |  Design  Sheet  ]  Curves  | 


Data:  |  FEA  Input  Data 


Name 

Value 

Units 

|  Description 

1 

Armature  Parallel  Branches 

1 

2 

Equivalent  Model  Depth 

135 

mm 

3 

Equivalent  Stator  Stacking  Factor 

0.95 

i 

Equivalent  Rotor  Stacking  Factor 

0.95 

5 

Equivalent  Br 

1.23 

tesla 

6 

Equivalent  He 

B90000 

A_per_meter 

7 

Region  Depth 

172.G37 

mm 

3 

Unit  Fractions 

2 

Figure  A-13  -  Solution  Data  Sheet  Showing  “FEA  Input  Data” 

The  Design  Sheet  tab  has  a  table  of  all  the  analysis  results.  The  Curves  tab  has  several  plots  to 
look  at. 
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The  D  and  Q  inductances  are  listed  as  “D-Axis  Inductance  Ll+Lad(H)”  and  “Q-Axis  Inductance 
Ll+Laq(H)”,  respectively,  and  they  are  found  under  “STEADY  STATE  PARAMETERS”. 

The  phase  resistance  is  listed  as  “Armature  Phase  Resistance  R1  (ohm)”.  It  is  also  found  under 
“STEADY  STATE  PARAMETERS”. 

If  the  number  of  conductors  per  slot  was  not  specified  (0  was  entered),  RMxprt  will 
automatically  determine  a  number.  This  number  can  be  found  under  “STATOR  DATA”  and  is 
listed  as  “Number  of  Conductors  per  Slot”.  So,  for  a  double  layer  concentric  winding 
arrangement,  half  of  those  conductors  will  belong  to  one  winding  and  half  to  another  in  the  same 
slot. 

Under  “STATOR  DATA”,  you  might  want  to  confirm  that  the  slots  were  not  over-filled: 

Limited  Slot  Fill  Factor  (%) :  90 

Stator  Slot  Fill  Factor  (%) :  38.143 

Using  RMxprt,  the  design  can  be  optimized  by  adjusting  various  inputs  and  quickly  seeing  what 
its  predicted  performance  would  be. 

Once  the  design  of  the  motor  has  been  finalized,  it  is  ready  to  convert  to  a  finite  element  method 
(FEM)  model.  The  following  section  explains  step-by-step  how  this  can  be  done  with  the  goal  of 
getting  three  of  the  five  parameters:  PM  flux  linkage,  direct  inductance,  and  quadrature 
inductance. 

A3. 3  Model  the  motor  with  finite  element  method  (FEM) 

A3. 3.1  Export  RMxprt  Model  to  Maxwell  2D 

A3. 3. 1.1  After  opening  and  selecting  the  RMxprt  design,  export  it  to  Maxwell  2D 


RMxprt  Tools  Window  Help 

1  List... 

Validation  Check.,. 

0  Analyze  All 
]g)  Edit  Notes... 

tp  #  la  1  ie  Iff*? 

Phase 

Turns 

In  Slot 

0i 

A 

1 

IT 

2B 

B 

1 

2T 

3B 

Machine  Editor 

-B 

1 

3T 

4B 

Design  Settings... 

■C 

1 

4T 

5B 

Analysis  Setup 

P  Add  Solution  Setup... 

List... 

Opti metrics  Analysis  ► 

Results  ► 

Export... 

Design  Properties... 

Create  Maxwell  Design... 

Design  Datasets... 

View  2D  Geometry 

Figure  A-14  -  Creating  a  Maxwell  Design  from  the  RMxprt  Model 
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A3. 3. 1.2  When  the  dialog  opens,  choose  2D 


Figure  A-15  -  Choosing  to  Create  a  2D  Maxwell  Design 

A3. 3. 1.3  To  see  what  the  new  design  looks  like,  double  click  on  the  design  name  under 
Project  Manager 

i-Jp  S12P10* 


Maxwell2D Design  1  (Transient,  XY'f 


|  i-^  Model 

i-O  MotionSetupl 

. M  Moving  1 

i  a-tf  Boundaries 
Masterl 
!■  Slavel 
L&  VectorPotentiall 
|  B-O  Excitations 
j  i-U|U  PhaseA 
|  PhAD 

i  td  Phil 

Figure  A-16  -  Project  Manager  Tree  for  Maxwell  2D  Design 
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A3. 3. 1.4  Set  solution  type  to  Magnetostatic 

1  -  JU  Modeler  -  [Uanaher  -  Maxwell^UUesignl  -  Modeler] 


Marvell  2D  Tools  Window  Help 


Solution  Type... 


IS  List... 

Validation  Check., 
|]  Analyze  All 
][il  Edit  Notes... 


■lodel 


3D  Model  Editor 

Design  Settings... 

Translate  Material  Database 


Model 

Boundaries 

Excitations 

Parameters 

Mesh  Operations 

Analysis  Setup 

Optimetrics  Analysis 

Fields 

Results 


Create  3D  Design. 


Export  Equivalent  Circuit 


Design  Properties... 
Design  Datasets... 


""U  r  u  IVIUVL 

■'□5  Rotate 


O  'i 


Figure  A-17  -  Selecting  “Solution  Type...”  from  “Maxwell  2D’ 
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Menu 


A3. 3. 2  Make  motor  whole 


A3. 3. 2.1  Select  the  user-defined  outer  vacuum.  Make  sure  it  is  under  OuterRegion 


l+j-L-Q  OUULrdLL 

vacuum 
[=]■■■□  Band 

| . &  CreateUserDefinedPart 

Subtract 
□■■■□  InnerRegion 

| . &  CreateUserDefinedPart 

□■■■B  Subtract 
□  OuterRegion 


. &  3 


••[Jl  CloneTo 
Shaft 

C  U  p  rfl  H  P  a  i+ 


Figure  A-18  -  Selecting  the  User-Defined  out  Vacuum 


A3. 3. 2.2  Change  the  fractions  to  1.  This  will  make  the  whole  motor  model  show 


Properties 


Name 

Value 

Unit 

Ev 

DLL  Locat... 

syslib 

DLL  Version 

12.0 

DiaGap 

43.1 

mm 

43.1mn 

DiaYoke 

77.1 

mm 

77.1mn 

Length 

0 

mm 

0mm 

SegAngle 

3 

deg 

3deg 

Fractions 

D 

1 

InfoCone 

m 

m 

Figure  A-19  -  Model  Fractions  Highlighted  in  the  Properties  Box 


Or  you  can  do  this  by  another  way.  Click  Maxwell  2D  on  top  of  the  menu  and  click  Design 
Properties. 
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A3. 3. 3  Add  Currents  and  Polarities 
A3. 3. 3.1  Select  the  coils  and  magnets 


Figure  A-20  -  Full  Model  with  all  Windings  and  Magnetic  Poles  Selected 
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A3. 3. 3. 2  Separate  the  bodies.  This  makes  it  so  you  can  modify  their  values  separately 


avell2DDesignl  -  3D  Modeler  -  [Danaher  -  Maxwell2DDesignl  -  Modeler] 
aw  [  Modeler  ]  Maxwell  2D  Tools  Window  Help 

g  m  p  m  I  ^ 


]  □+ 

Import  From  Clipboard 


- [s  Group  Objects  By  Material 


,n^  ^  jFf  Assign  Material... 


Import... 

Export... 


nit 


Movement  Mode 
Snap  Mode... 
New  Object  Type 


Coordinate  System 
List 


Edge 

Surface 


Boolean 


Units... 

Measure 


Generate  History 

Delete  Last  Operation 
Purge  History 


Model  Analysis 


^  |  |  Model 


3 


Ar 


^  LLnite 
CJ  Subtract... 
e  Intersect 
□□  Split- 


Separate 


rt 

Ar 


3-%  Unite 


Figure  A-21  -  Separating  the  Windings  and  Magnets  so  that  they  are  no  Longer  Linked 
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A3. 3. 3. 3  Add  variables  for  phase  currents 


Project  Manager 


3-Jp  Danaher" 


a-  m  Maxwell2D 


Dfisinnl  fMannfth-ia^ir:  XtT 


+ 


6?  Model 
^  Bound. 
®  Excitat 
^  Param* 
Mesh 

JF  Analysi 
Lj^l  Optime 
■|S|  Result: 


Project  f 


rrr 


Properties 


Name  Value 


Copy 

Paste 


X  D  * 


Rename 

Delete 


H 


Solution  Type... 
List- 

Validation  Check... 

Analyse  All 
Edit  Notes... 


3D  Model  Editor 

Design  Settings... 

Translate  Material  Database  , 


Create  3D  Design., 


Export  Equivalent  Circuit 


Design  P^p 


nerties., 


Design  Datasets... 


:'""l 


IT 

Figure  A-22  -  Selecting  “Design  Properties...”  from  “Maxwell2DDesignl  (Magnetostatic 

XY)” 


A  recommended  set  of  variables  is  shown  here: 
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C  Optimization 


C  Tuning 


C  Sensitivi 


Local  Variables 


(*  Value 


Name 

Value 

Unit 

1 

1 

A 

1A 

th 

0 

0 

thRad 

{th*pi/1S0} 

0 

phiShrft 

0 

0 

phi  Rad 

ptiShftWISO 

0 

Hi  120 

pi-2/3 

209439. 

iaa 

Tcos#hRad  +  phi  Rad} 

1A 

ibb 

Tcos^hRad  +  phi  Rad  -  th  1 20} 

-0.5A 

icc 

rcosfthRad  +  phi  Rad  +  Uhl  20} 

-0.5A 

Nc 

m 

36 

iaa Nc 

iaa'Nc 

36 

ibb Nc 

ibb'Nc 

43 

icc Nc 

icc'Nc 

43 

P 

10 

10 

thM 

thRad-2/P 

0 

Figure  A-23  -  Recommended  List  of  Design  Variables 


i 

current  amplitude 

th 

electrical  angle  of  rotor  [degrees] 

thRad 

electrical  angle  of  rotor  [rads] 

phiShift 

torque  angle  [degrees]  (atan(iq/id)) 

phiRad 

torque  angle  [rads] 

th!20 

120  degrees  converted  to  radians.  P 

ACB  winding. 

iaa 

phase  A  current 

ibb 

phase  B  current 

icc 

phase  C  current 

Nc 

turns  per  coil 

iaa_Nc 

total  current  thru  a  phase  A  coil 

ibb  Nc 

total  current  thru  a  phase  B  coil 

icc_Nc 

total  current  thru  a  phase  C  coil 

P 

number  of  poles 

thM 

mechanical  angle  of  rotor  [r 

Note  that  when  you  create  these  variables,  only  the  current  amplitude  /  must  have  a  specified 
unit.  You  do  this  by  typing  the  unit  after  its  nominal  value: 
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- ■  '  ■ 

Add  Property 


Name  f[ 


(•  Variable 

C  Value  C  Men 


Value 


14 


Enter  a  number  with  units,  variable,  or  expression  into  the  Value  field.  Referer 
$C1 , 2“cos($x). 


Figure  A-24  -  Defining  the  Nominal  Value  of /with  the  Unit  A  (Ampere) 

A3. 3. 3. 4  Add  excitations.  Consider  that  it  might  be  helpful  to  refer  to  the  RMxprt  winding 
layout  before  creating  the  coil  excitations.  It  could  also  be  helpful  to  first  color 
code  your  coils  to  guide  you  in  assigning  currents. 


Figure  A-25  -  An  Image  from  RMxprt  of  the  Winding  Layout  for  the  Motor  (Left)  and  an 
Image  from  Maxwell  2D  of  the  Color-coded  Winding  Layout  (Right).  The  Colors  Must  be 

added  by  the  User 


When  you  are  ready  to  create  excitations,  select  one  coil  at  a  time  and  assign  currents.  The 
reason  for  assigning  current  to  one  coil  at  a  time  is  for  control  over  naming. 
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Make  sure  the  value  is  the  phase  current  times  the  number  of  turns  per  slot,  Nc.  Note,  if  the  coil 
is  a  return  path,  do  not  put  in  a  negative  value,  but  simply  select  Negative  for  Ref.  Direction. 
Also,  note  that  the  excitations  on  the  coils  should  be  symmetrical. 


Name;: 


|  PhA 


Parameters 

Value: 


iaa  Nc 


[a  3 


Ref,  Direction:  f*  Positive  C  Negative 


Figure  A-27  -  Setting  the  Current  Excitation  for  a  Winding 

For  a  concentric  winding,  it  is  worth  naming  a  coil  and  its  return  path  by  very  similar  names  (e.g. 
PhA  l  and  PhARel).  This  is  why  it  is  best  to  assign  currents  to  one  coil  at  a  time. 

Assigning  excitations  is  a  very  delicate  process.  It  is  very  important  to  double  check  your 
assignments.  It  would  help  to  check  that  the  phase  names  match  the  color  coding  and  are  in  the 
right  order;  then  check  that  the  current  values  and  Ref.  Directions  are  appropriate  for  the 
excitation  names;  then  check  that  each  coil  excitation  is  properly  named  to  match  its  return  path 
coil. 
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A3. 3. 4  Set  magnetic  polarities 


Assign  Magnetic  material  with  polarity.  Unfortunately,  when  the  motor  model  was  made  whole, 
Ansys  made  the  other  side  a  mirror  image,  which  depending  on  the  number  of  poles  means  the 
magnets  on  the  lower  half  will  have  the  wrong  polarity. 


!--|P  Optimetrics 
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□■■■□  Mag_l 
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E-IB  Subtract 
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□■■■n  Band 
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Figure  A-28  -  Selecting  the  Polarity  of  Magnets 


A3. 3. 5  Assign  vector  potential 


A3. 3. 5.1  Select  Edges  of  the  outer  edge  of  motor 


Figure  A-29  -  Setting  the  Selection  Option  to  Edges 
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A3. 3. 5.2  Assign  Vector  Potential  to  outer  edge.  The  value  should  be  0  weber/m 

~  *  * 


Project  Manager 


El  JO— I  Danaher* 

□  SI  Maxwell2DDesign1  (Magnetostatic,  J 
&  Model 
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Vector  Pot 
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rrr 


Project  r 
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:eFro 


an 

l±H 


□  °  Rotate 
NdFe30_N 
NdFe30  S 


Figure  A-30  -  Selecting  “Vector  Potential...”  from  “Assign”  from  “Boundaries’ 

A3. 3. 6  Analyze 

A3. 3. 6.1  Add  solution  setup.  You  can  use  all  the  default  settings. 


Project  Manager 


3  -ftJ  Danaher* 

H  □!  Maxwell2DDesign1  (Magnetostatic^ 

j . &  Model 

Boundaries 
S-®  Excitations 
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Mesh  Operations 
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Project  j~ 
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List... 
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Figure  A-31  -  Selecting  “Add  Solution  Setup...”  from  “Analysis’ 
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A3. 3. 6.2  Assign  Matrix.  A  matrix  is  used  to  define  how  the  various  coils  are  connected  and 
what  makes  a  loop. 
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Figure  A-32  -  Selecting  “Matrix...”  from  “Assign”  from  “Parameters” 

A3. 3. 6. 3  Include  currents,  but  not  all.  In  the  Matrix,  specify  the  return  path  for  each 

conductor.  If  one  conductor  is  selected  to  be  the  return  of  the  other  conductor, 
don’t  include  it  in  the  matrix.  For  example,  don’t  include  PhaseARe  l  in  the 
matrix.  Then  for  the  return  coil,  make  its  return  path  be  “infinite”. 
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Figure  A-33  -  Selecting  all  but  Return  Windings 

A3. 3. 6.4  Group  currents  into  phases.  Make  sure  to  specify  the  number  of  turns.  Also,  it  is  a 
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good  idea  to  name  your  groups  (e.g.  PhA,  PhB,  PhC). 
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Figure  A-34  -  Grouping  Windings  into  Phases 


A3. 3. 6. 5  Validate  the  Model 


[Maxwell  2D  Tools  Window  Help 

Solution  Type... 

IS  List... 

Validation  Check... 

[|]  AnalyzeAll 
][i|  Edit  Notes... 

3D  Model  Editor 

Figure  A-35  -  Selecting  “Validation  Check...”  from  “Maxwell  2D”  Menu 

A3. 3. 6. 6  Analyze  the  Analysis  Setup 

ij  -ffl  Mesh  Operations 


Analysis 
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Revert  to  Initial  Mesh 

Figure  A-36  -  Selecting  “Analyze”  from  “Setupl” 
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A3. 3. 6. 7  Look  at  Solution  Data 
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Figure  A-37  -  Selecting  “Solution  Data...”  from  “Results’ 


A3. 3. 6. 8  Look  at  inductances  and  fluxes.  Note,  until  you  select  the  “PostProcessed”  check 
box,  the  inductances  are  for  each  coil.  Keep  in  mind  that  these  inductances  are 
actually  inductances  per  unit  meter.  Therefore,  these  values  must  be  multiplied  by 
the  length  of  the  motor  to  get  the  true  inductances  of  the  motor. 


Figure  A-38  -  Post-processed  Inductances  for  Phases  A,  B,  and  C 
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There  is  a  Matlab  function  made  for  this  analysis  called  “dqO.m”.  It  will  automatically  convert 
arrays  in  ABC  reference  frame  to  DQO  reference  frame.  In  the  case  of  the  example  in  the  figure 
above,  the  motor  length  is  0.058  m.  Here  you  can  see  how  the  first  argument  to  dqO  is  the  ABC 
matrix  of  mutual  and  self-inductances,  the  second  argument  is  the  angle  of  the  rotor  converted  to 
radians,  and  the  third  parameter  is  the  indicator  that  this  is  an  ABC  wound  motor  (as  opposed  to 
ACB  wound).  Notice  that  the  result  is  multiplied  by  0.058  to  account  for  the  length  of  the  motor. 
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17 .273263333232321 
0.000000424090499 


-0 .263350963304730 
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Figure  A-39  -  Matlab  Command  Window  Showing  the  Results  of  Using  dqO  on  the  ABC 

Phase  Inductances 


So,  the  top  left  value  in  the  result  is  Ld ,  the  middle  value  is  Lq ,  and  the  bottom  right  value  is  L0 . 
Just  because  L0  is  not  zero  does  not  mean  that  the  motor  is  not  balanced.  Also,  note  that  the 
results  from  Maxwell  2D  were  in  mH  so  the  results  from  Matlab  are  also  in  mH.  This  indeed 
gives  valid  inductance  values;  however,  inductance  is  a  function  of  currents.  So,  a  following 
section  will  discuss  how  to  get  the  inductances  as  a  function  of  currents. 

To  get  the  PM  flux  linkage,  the  test  must  be  run  with  current  set  to  zero  the  whole  time.  Maxwell 
will  calculate  the  flux  linkages  into  the  windings.  We  want  these  values  to  be  the  result  of  the 
permanent  magnets  only,  not  from  any  flux  set  up  by  neighbor  windings.  Furthermore,  the  flux 
linkage  from  the  permanent  magnet  should  be  independent  of  current. 
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S3  Solutions:  Danaher  -  Design2B 


i  „H5ib*5*d' 


Simulation: 


|Setup2 


^  |  |  LastAdaptive 


Design  Variation:  |_Nc='[r  Nc='8G'  P='10'  phiRad='0'  phiShift='0'  th='0'  thl  20=2. 0943951 023932'  thM='0'  thRad='0' 


Figure  A-40  -  Flux  Linkages  for  Phases  A,  B,  and  C 

In  much  the  same  was  as  with  inductances,  the  flux  linkages  can  be  converted  to  DQO  reference 
frame  with  Matlab. 


Command  Window 


i+-  n  ?  x 


»  dqO( [2.5788,  -1.2562 
an  s  = 


|  ^  Editor  -  edract lambdapm.m 
,  -1.2823],  0 . 2*pi/13Q, 


Figures  -  Figur 


true ) *0 . 053 


0.143793410933637  0.000354609067273  0.000779133333333 


h  »  I 

Figure  A-41  -  Matlab  Command  Window  Showing  the  Result  of  Using  dqO  on  the  ABC 

Phase  Flux  Linkages 

Notice  that  again,  the  values  must  be  multiplied  by  the  length  of  the  motor  (in  this  case  0.058  m). 
The  far  left  value  in  the  result  is  the  flux  linkage  in  the  D  axis.  This  is  the  flux  linkage  of  the 
permanent  magnet.  Now,  these  flux  linkages  are  not  inherently  particular  to  the  permanent 
magnet.  Only  because  the  currents  are  zero,  the  only  flux  through  the  stator  windings  is  due  to 
the  permanent  magnets.  You  will  notice  that  the  units  in  Maxwell  2D  are  probably  Wb.  So,  the 
result  from  Matlab  is  also  in  Wb.  Unfortunately,  there  does  not  seem  to  be  any  way  to  confirm 
this  value  with  RMxprt. 

The  next  section  shows  how  to  get  the  nonlinear  inductances  as  functions  of  currents. 
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A3. 4  Extract  the  reduced-order  parameters 

The  purpose  of  using  parametrics  is  to  try  many  combinations  of  different  settings  such  as 
different  values  of  current,  phase  angle,  rotor  angle,  magnet  thickness,  etc.  The  following  steps 
show  exactly  how  these  parametrics  are  setup  and  then  simulated. 

A3. 4.1  Simulate 

A3. 4. 1 . 1  Add  parametrics 


Js]  r~ 

SF1 

|  0O  C 

i _ 

■S  Result 

@  Paste 

|  i  C 

Reid  £ 

Add 

v_1  k. 

v=a  Parametric...  ^ 

RMxprtDe 
|  Definitions 

Analyze  ► 

Juning... 

Vi  ew  An  a  lys  is  Res  u  It. . . 

X  Optimization... 

Sensitivity...  ^ 

[||  Statistical...  t 

■'  in  :  r  sA 

Figure  A-42  -  Selecting  “Parametric...”  from  “Add”  from  “Optimetrics” 
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A3.4.1.2  Click  Add...,  select  one  of  the  variables  from  the  list,  specify  the  pattern  type  (e.g. 
Linear  step),  specify  the  Start,  Stop,  and  Step  values,  and  then  the  “Add  »” 
button.  Make  sure  to  press  this  button  or  else  nothing  will  happen.  To  get  a  good 
sampling  of  the  inductances,  make  sure  to  vary  current  magnitude  from  near  zero 
to  some  positive  number,  torque  angle  from  0  to  pi,  and  rotor  angle  from  0  to 
some  fraction  of  2*pi.  The  user  will  need  to  test  in  order  to  find  the  ranges  that  are 
best  for  the  particular  motor. 


Setup  Sweep  Analysis 


£3 


Sweep  Definitions  |  Table  )  General  ]  Calculations  ]  Options  ] 


Sync  tt  Var  f 


Operation 


Add/Edit  Sweep 


Variable  \\ 

C  Single  value 
'■*  Linear  step 
r  Linear  count 
f  Decade  count 
r  Octave  count 
r  Exponential  count 


~3 


Start: 


Step:  [T 


F — 

|A — 3 

[20 

F  3 

r 

|a“3 

Add  »  | 

Update  >>  | 
Delete 


Variable 


Description 


Linear  Step  from  1A  to  2QA,  s 


OK 


Cancel 


OK 


Add.. 


Edit... 


Delete  | 


Cancel 


Figure  A-43  -  Setting  up  a  Sweep  Analysis 

A3. 4. 1.3  Now  analyze  from  the  parametric  setup.  This  could  take  quite  a  while  depending 
on  how  many  parametric  combinations  you  set  up. 


□■■■ijg]  Optimetrics 

! . _ 

j . |§]  Results 

Field  Overlays 

^  Rf^xprt  Design  1  (Brushlesi 
n  Definitions 

Hi  Copy 
@  Paste 

Rename 

X  Delete 

rrr 

J 

Properties... 

Analyze 

s 

Vi  ew  An  a  lysi  s  Res  u  It. . . 

Figure  A-44  -  Selecting  “Analyze”  from  “Parametric  Setup  1” 
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A3. 4. 1.4  For  numerical  analysis  on  the  results,  create  a  data  table. 


Fielc 
RMxprtE 
^  Definitioi 


Paste 


Create  Magnetostatic  Report 


Create  Fields  Report 


Delete  All  Reports 


[j-l=]  C 

a-n  C 


Rectangular  Plot 


Data  Table 


r 


3D  Rectangular  Plot 

1 


OT 


Figure  A-45  -  Selecting  “Data  Table”  from  “Create  Magnetostatic  Report”  from  “Results’ 


A3.4.1.5  For  an  inductance  report,  select  all  the  mutual  and  self-inductances  in  the  list. 


Trace  |  Families  | 


Solution:  | Setup  j  .  LastAdaptive 

Primary  Sweep:  |l  |  All 

Parameter:  |None 

~3 

X:  W  Default  fi 

Y:  [Matrix !.L(PhA,PhA);  Matrixl.L(PhA,PhB);  Matrix l.L(PhAfPhC);  Matrix !.L(PhB,PhA);  Matrix !.L(PhB,PhB);  Matrix  ^ 

Function: 

abs 

acos 

acosh 

ang_deg 

ang_rad 

asin 

asinh 

atan 

atanh 


Figure  A-46  -  Selecting  all  Mutual  and  Self-inductances 


A3. 4. 1.6  For  a  PM  flux  linkage  report,  make  sure  to  select  MagFlux  and  all  the  phases  in 
the  list. 


Context 


Solution: 

Parameter: 


j  Setup  1  :  LastAdaptive 


None 


31 


Trace  |  Families  | 

Primary  Sweep:  F - nr  All 

X:  |7  Default  |1 

Y:  |Matrixl.MagFlux(PhA);  Matrix !.MagFlux(PhB);  Matrix  !.MagFlux(PhC) 


Category: 

Variables 
|  Output  Variables 
L 

ISBSHTffl 

MagFluxNom 

Design 


Figure  A-47  -  Selecting  all  Magnetic  Flux  Linkages 
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A3. 4. 1.7  Once  the  reports  are  generated,  transpose  the  data  so  that  the  long  labels  are  in  a 
column  on  the  left. 


•••■ 


Data  Table  1 


■  g  Matral.UPhA.PhA) 

■  g  Matral.UPhA.PhB) 

•  g  Matral.UPhA.PhC) 

•  g  Matral.UPhB.PhA) 
g  Matral.UPhB.PhB) 

■  g  Matral.UPhB.PhC) 

■  g  Matral.UPhC.PhA) 

•  g  Matral.UPhC.PhB) 

•  g  Matral.UPhC.PhC) 


■in  Data  Table  2 

g  Matrix  1 .  Mag  FluxfPh  A) 
■g  Matra  1.  Mag  FlLK(PhB) 
•g  Matrix  1.  Mag  RuxlPhC} 


n i  j  r* i. 


2 

2.000000 

3 

3.000000 

i 

a  nnnnnn 

Modify  Report ... 

6 

Update  Report 

7 

Font ... 

S 

9 

Print 

10 

Export  Data... 

11 

Import  Data... 

12 

13 

Copy 

U 

Transpose 

IS 

Reset  Order 

Ifi 

Accumulate 

Figure  A-48  -  Transposing  the  Data  Table 


A3.4.1.8  Copy  the  data  by  right  clicking  in  the  table 


i.oooot 


|hift='0l  th='0' 

Modify  Report ... 

Update  Report 

22.0930* 

|hift='0'  th='0' 

-2.7073; 

iiShift='0'  th='0' 

Font ... 

-2.6696; 

|iShift='0'  th='0' 

Print 

-2.7073; 

iiShift='0'  th='0' 

22.1188; 

|iShift='0'  th='0' 

Export  Data... 

-2.6788* 

iiShift='0'  th='0' 

Import  Data... 

-2.6696; 

jiS  hift='0'  th='0' 

Copy 

-2.6788* 

jjShift='0'  th='0' 

Transpose 

22.2142* 

Figure  A-49  -  Copying  the  Data 
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A3.4.1.9  This  can  be  then  pasted  into  Excel.  Save  this  as  a  csv  file 


&  Cut 
Copy 

if  Format  Painter 
Clipboard  ^ 


Font: 


A1 


fir  I [A] 


A  B 
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E 
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i 
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3 

4 

2 
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21 

3 

Matrix  l.L(  -2.70734 

-2.67386 
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-2.67358 

-2.6: 

4 

Matrix l.L(  -2.66963 

-2.64767 

-2.65418 

-2.65484 

-2.E 

5 

Matrix l.L(  -2.70734 

-2.67386 

-2.67385 

-2.67358 

-2.6: 

6 

Matrix  l.L(  22.11884 

21.92841 

21.94362 

21.95  782 

21.9: 

7 

Matrix  l.L(  -2.67886 

-2.64693 

-2,64718 

-2.64736 

-2.6^ 

8 

Matrix l.L(  -2.66963 

-2.64767 

-2.65418 

-2.65484 

-2.E 

9 

Matrixl. L(  -2.67886 

-2.64693 

-2.64718 

-2.64736 

-2.6^ 

10 

Matrix  l.Lf  22.2142 

21.9997 

22.00232 

22.02007 

22.0^ 

Figure  A-50  -  Pasting  the  Data  into  Excel 

Now,  the  data  is  safely  in  a  spreadsheet,  but  notice  that  depending  on  how  many  simultaneous 
parametric  variables  you  have,  the  distinction  of  what  parameters  had  what  values  is  all  together 
in  one  column  on  the  left. 

To  extract  the  data  into  a  more  usable  format,  we  need  to  run  a  Matlab  script.  This  process  is 
explained  in  the  following  section. 

A3. 4.2  Extract  the  data 

Here  are  two  typical  lines  in  the  first  column  of  one  of  the  csv  files  we  might  generate: 

Matrixl .L (PhA, PhA)  [mH]  -  I='64A'  Nc='86'  P='10' 
phiShift=' 180.1 '  thl20='2. 0943951023932'  -  Setup2  : 

Last Adaptive 

Matrixl .L (PhA, PhA)  [mH]  -  1= ' 45 . 254833995939A'  Nc='86' 

P='10'  phiShift=' 180 . 1 '  thl20= ' 2 . 0943951023932 '  -  Setup2  : 
La st Adaptive 

These  text  strings  contain  the  labels  of  the  particular  test  that  was  run,  but  we  need  to  parse  the 
data  into  a  multidimensional  array  that  can  be  indexed  based  on  each  of  the  relevant  labels  (e.g. 
phase  pair  “PhA, PhA”,  current  value  “I=’45.25”\  phase  shift  “phiShift=’ 180.1’”).  This 
multidimensional  array  will  be  much  easier  to  analyze  than  the  csv  file.  To  parse  the  data  into 
such  an  array  we  need  the  provided  Matlab  script  csvlabeled2num.m.  This  script  will  need  to 
know  what  labels  to  look  for.  So,  first  we  need  to  define  the  labels  as  cell  arrays. 

phiCell  =  { 'phiShift=' '0.1' ,  'phiShift= ' ' 10 . 1 ' , 

' phiShift= ''20.1',...  ' phiShif t= ' ' 30 . 1 ' ,  'phiShift= ' ' 40 . 1 ' , 
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' phiShif t= ''50.1',...  ' phiShif t= ' ' 60 . 1 ' ,  'phiShift= ' ' 70 . 1 ' , 

' phiShift= '  '80.1',...  ' phiShif t= '' 90 . 1 ' ,  ' phiShift= '  ' 100 . 1 '  , 

' phiShif t= ' '110.1',...  ' phiShif t= ' '120.1', 

' phiShif t= ' '130.1',  'phiShift=' '140.1',... 

' phiShif t= ' '150.1',  ' phiShif t= ' '160.1', 

' phiShif t= ' ' 170 . 1 ' , . . .  'phiShift=' ' 180 . 1 ' } ; 
currentCell  =  {'I=''1A',  'I=''1.4',  'I=''2A',  ' 1= ' '2.8', 

' !=' ' 4A' ,  ... 


' 1= '  '5.6',  ' 1= '  ' 8A ' ,  'I=''ll',  ' 1= '  '16',  '  1=  '  '  22  '  , 

' 1= ' ' 32 ' ,  ...  ' 1= ' ' 45 ' ,  ' 1= ''64'}; 

abcCell  =  { ' PhA, PhA ' , ' PhA, PhB ' , ' PhA, PhC ' , . . . 

' PhB, PhA' , ' PhB, PhB' , ' PhB, PhC' , . . . 

' PhC , PhA ' , ' PhC , PhB ' , ' PhC , PhC ' } ; 


Notice  the  double  quotes  ‘  ’  used  to  create  a  single  quote  ‘  in  the  actual  string. 

Then,  these  cell  arrays  need  to  be  used  in  the  csvlabeled2num.m  script. 

theArray  =  csvlabeled2num ( ' LL  vs  phi,  i,  abc,  theta . csv ',.. . 
{phiCell,  currentCell,  abcCell }) *motorLength; 

The  script  takes  two  parameters:  the  name  of  the  csv  file  and  a  cell  array  of  the  cell  arrays. 

Notice  that  the  cell  arrays  must  be  in  the  order  that  you  want  the  data  organized  in  dimensions  in 
the  array.  So,  in  the  above  example,  the  inductance  values  will  vary  with  phi  along  the  first 
dimension,  with  current  along  the  second  dimension,  and  with  ABC  pairing  along  the  third 
dimension.  The  last  dimension  is  always  the  actual  second  dimension  in  the  csv  file.  So,  in  the 
above  example,  each  column  in  the  csv  file  was  for  a  different  value  of  theta.  Thus,  the  fourth 
dimension  varies  theta  in  the  array. 

This  analysis  can  quickly  become  rather  complicated.  So,  to  simplify  the  problem,  use  the 
Inductances  package.  To  extract  the  inductances,  run  extract.m.  The  inputs  are  set  in  the  section 
labeled  "Get  the  array  of  inductances  from  the  csv  file."  There  are  two  output  files:  L(phi,i).mat 
and  L(phi,i,theta).mat.  The  inductances  from  a  typical  motor  would  look  something  like  these: 
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Figure  A-51  -  Direct  inductance  (left)  and  quadrature  inductance  (right),  both  as  functions 

of  direct  and  quadrature  currents 

To  smooth  the  inductances  and  calculate  their  derivatives  with  respect  to  id  and  iq,  run 
smooth2d.m.  Its  input  is  L(phi,i).mat.  It  smooths  the  inductances  and  calculates  the  derivatives 
with  respect  to  id  and  iq.  Its  output  file  is  L(phi,i),  DL(id,iq)  smooth.mat.  After  smoothing,  the 
inductances  would  look  more  like  these: 


Figure  A-52  -  Direct  (Left)  and  Quadrature  (Right)  Currents  Smoothed 


These  tables  of  inductances  then  need  to  be  placed  into  the  motor  module  in  the  EMA  package, 
in  the  Modules  folder. 

In  addition  to  developing  the  inductances,  the  FEM  model  that  was  built  for  Maxwell  2D  can 
also  be  ported  into  other  ANSYS  software  to  do  thermal  analysis.  Here  short  FEM  simulations  of 
thermal  response  can  be  run  to  collect  data  for  developing  lumped  element  parameters  of  thermal 
resistance  and  thermal  capacitance. 
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A4  Conclusion 


While  many  of  the  necessary  parameters  for  EMA  simulation  can  be  found  through  non¬ 
destructive  experiments,  some  parameters  are  best  determined  by  using  a  closely  matched  model. 
The  main  parameters  that  need  to  be  found  by  modeling  are 

•  PM  flux  linkage 

•  Direct  inductance 

•  Quadrature  inductance 

•  Thermal  capacitance 

•  Thermal  resistance 

The  process  to  find  these  requires  great  care  and  precision,  but  they  can  be  found  without 
destroying  the  EMA. 
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LIST  OF  ACRONYMS,  ABBREVIATIONS,  AND  SYMBOLS 


ACRONYM 

DESCRIPTION 

AFRL 

Air  Force  Research  Laboratory 

EMA 

Electromechanical  Actuator 

EU 

Electronic  Unit 

FEM 

Finite  Element  Method 

FOC 

Field-oriented  Control 

HPEAS 

High  Performance  Electric  Actuation  Systems 

IGBT 

Insulated-gate  Bipolar  Transistor 

INVENT 

Integrated  Vehicle  and  Energy  Technology 

NL-LEM 

Nonlinear  Lumped-element  Model 

PCM 

Phase  Change  Materials 

PM 

Permanent  Magnet 

PMSM 

Permanent  Magnet  Synchronous  Machine 

PWM 

Pulse-width  Modulation 

UCF 

University  of  Central  Florida 
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